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Dispersive representations of the tttv scattering amplitudes and pion form factors, valid at two-loop 
■ accuracy in the low-energy expansion, are constructed in the presence of isospin-breaking effects 

induced by the difference between the charged and neutral pion masses. Analytical expressions for 
• the corresponding phases of the scalar and vector pion form factors are computed. It is shown 

[jL^ , that each of these phases consists of the sum of a "universal" part and a form-factor dependent 

' contribution. The first one is entirely determined in terms of the irn scattering amplitudes alone, 

£Nj ■ and reduces to the phase satisfying Watson's theorem in the isospin limit. The second one can 

be sizeable, although it vanishes in the same limit. The dependence of these isospin corrections 
with respect to the parameters of the subthreshold expansion of the nir amplitude is studied, and 
O-i an equivalent representation in terms of the S-wave scattering lengths is also briefly presented 

and discussed. In addition, partially analytical expressions for the two-loop form factors and nir 
D , scattering amplitudes in the presence of isospin breaking are provided. 



I. INTRODUCTION 



In recent years, our knowledge of low-energy pion-pion scattering has improved in a very significant way 
and in several respects. Firstly, the high precision K^ 4 experiments performed at the BNL AGS by the E865 
experiment [H, 0] and, more recently, by the NA48/2 collaboration 0,0] at the CERN SPS, have provided very 
accurate determinations of the difference Sq — 5\ of the pion-pion phase shifts in the S and P waves in the 
energy range between threshold and the kaon mass. Next, one should mention the measurement of the invariant 
mass distribution in K 7r ± 7r°7r° decays [a, |6(, that gives information on the S- wave 7T7t scattering lengths 
(see also jsTtloll). Finally, forthcoming analyses of the data collected by the NA48/2 experiment on the K^ 4 
^ . decay channel into a pair of neutral pions (for preliminary reports, see [ill E3) 



decay channel into a pair of neutral pions (for preliminary reports, see [ll|,ll2j), or on the K\ — > tt tt tt decay 



mode [13|, together with the measurement of the pionium lifetime by the DIRAC collaboration [14j, should 
provide additional information, and might sharpen the picture even more. In the meantime, the accuracy 
obtained on K14 decays from NA48/2 implies that these data clearly drive the current determination of the 
difference between the S and P phase shifts at low energies, and in particular of the two scattering lengths 
<Zq and <2q, for which very accurate predictions are available [15|. This provides a particularly stringent test of 



two-flavour chiral perturbation theory [lq . and its underlying assumptions [T7I ] . 

In order to extract relevant information on low-energy pion-pion scattering from the above processes, it 
has become mandatory to take isospin violations into account. This is certainly quite easy to understand in the 
case of the — > ^tt ^ decay, where one exploits the presence of a unitarity cusp in the invariant tt°tt° mass 
distribution, which occurs only if the masses of the charged and neutral pions differ 0, ■ Perhaps somewhat 
more unexpectedly, isospin- violating corrections proved also of importance ,19] in the analysis of the K^ 4 data, 
in order to account for the high precision reached by the recent NA48/2 experiment, and to make comparison 
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with theory meaningful 2CJ, |2l| . Actually, once isospin corrections are applied also to the E865 data, there 



remains a disagreement with NA48/2 20(, whose origin seems to lie in the original analysis performed by the 
E865 collaboration (for details, see the errata quoted under refs. (HQ)- Anyway, the analysis of the full data set 
collected by NA48/2 has by now completely superseded the E865 results, and one should focus on the former 
to study pion-pion scattering from decays. 

In the present paper, we propose to address the issue of isospin- violating effects in low-energy pion-pion 
interactions using an approach based on a dispersive construction of the various tttt scattering amplitudes and 
pion form factors in the presence of isospin breaking. Ultimately, we wish to extend this program [22j to the 
Ke4 form factors analysed in the NA48/2 experiment. Before undertaking this enterprise and hitting the full 
complexity of this four-body decay, we want to demonstrate its feasibility and exhibit the general features of 
such a method by considering the somewhat simpler setting provided by the scalar and vector form factors of 
the pion. 

As far as the amplitude for elastic tttt scattering in the isospin limit is concerned, the general framework 
has been laid down in ref. [23j , and the explicit construction of the two- loop amplitude has subsequently been 
performed along these lines in detail in ref. 24|. Concerning the pion form factors, the corresponding dispersive 
representations in the framework of the chiral expansion have been studied in ref. in the isospin limit, but 
only one-loop expressions were given in analytical form. Full two-loop expressions of the vector form factors 



have been obtained by integrating the corresponding dispersive integrals in ref. |26| . In ref. |27j the two- loop 
expressions of the vector and scalar form-factors have also been obtained in the absence of isospin violation 
by the direct evaluation of Feynman graphs generated from the effective chiral Lagrangian at next-to-next-to- 
leading order. A similar calculation for the pion-pion scattering amplitude in the isospin limit has been achieved 
in ref. 28j. Finally, let us also mention that the reconstruction theorem for elastic tttt scattering in the isospin 



limit of ref. |23[ was extended by the authors of ref. [29( to the whole set of scattering amplitudes involving the 
mesons of the lightest pseudoscalar octet. Applications of this framework to the decay modes P — > tttttt, with 



P = K,rj 1 have also been considered 30, 31] 



These dispersive constructions generate subtraction polynomials with unspecified coefficients. The latter 
are in one-to-one correspondence with the appropriate combinations of low-energy constants and chiral loga- 
rithms that would be encountered in a calculation of the corresponding Feynman diagrams generated by the 
chiral lagrangian. In the case of the form factors, these coefficients may be identified with their slopes and 
curvatures. In the case of the tttt scattering amplitudes, they can be expressed in terms of the subthreshold 
parameters occurring in the expansions of these amplitudes as Taylor series around the center of the Mandel- 
stam triangle. This was the option considered in the isospin-symmetric case in ref. 24|. By no means, however, 
is this choice a necessity. It has, for instance, become customary to rather let the scattering lengths play a 
prominent role. They have a more direct physical interpretation than the subthreshold parameters, and are thus 
considered as more "experimentalist friendly" . We will therefore also provide expressions where the subtraction 
polynomials are given in terms of the two S-wave 1 = and 1 = 2 scattering lengths, a[j and aL in the isospin 
limit. In the isospin-symmetric situation, this provides an alternative to the choice made in ref. 24] ■ Of course, 
taking the expressions at two-loop order provided in the latter reference, one could convert the expressions 
for the tttt scattering amplitude given there to the one presented here in terms of the scattering lengths. The 
two formulations are equivalent, up to corrections that are of higher order. In the situation where isospin is 
broken, this allows us to discuss the size of the corresponding corrections to the phases of the form factors in 
terms of Oq and <2q. This second option is of course the most interesting in the present context, where these 
scattering lengths are the quantities one would eventually like to determine from the data. It is thus important 
that the corrections due to isospin breaking are not studied for a fixed a priori value for them. Indeed, given 
the precision reached by the latest experiments, one ought to perform a quantitative evaluation of the possible 
bias introduced if isospin corrections are evaluated for fixed values of these scattering lengths. This provides 
another motivation for the present work. 

Here we will mainly concentrate on the phases of the pion form factors. Full two-loop expressions for the 
scattering amplitudes and form factors themselves require the evaluation of dispersion integrals corresponding 
to specific topologies of two- loop three-point Feynman diagrams of the non- factorizing type ( "acnode" of "fish" 
diagrams, cf. fig. Explicit analytical expressions for them do not seem to be available in the literature in 
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the cases where several distinct masses are present. We therefore present only partially analytical expressions 
for the scattering amplitudes and form factors. Note that a similar situation arises in the evaluation of the 
SU(3) vector 32 1 and scalar [3^] form factors at two loops without isospin breaking, but where the difference 
between the pion and kaon masses has to be dealt with. These difficulties do not show up in the computation 
of the phases of the two-loop form factors, where the technically most demanding step is the computation of 
the projections of the one-loop amplitudes on the S and P partial waves, which can be done analytically. 

Coming now to the outline of this paper, our first purpose will be to extend the frameworks of refs. 



23, 25] to the situation where the difference between the masses of charged and neutral pions is taken into 



account. The general framework leading to two-loop representations for form factors and scattering amplitudes 
when isospin symmetry no longer holds is thus described in section [TTJ In section HTT1 we implement the program 
of constructing the corresponding form factors and scattering amplitudes at the one-loop level and provide 
explicit expressions for them. The second iteration, leading to two-loop representations of the form factors and 
the scattering amplitudes, is discussed in section |IV] The issue of isospin breaking in the phases of the two-loop 
form factors is addressed in section [V] Section I VII is devoted to the numerical analysis of the isospin-breaking 
contributions in the phases of the form factors. Finally, a summary of this study and our conclusions are to be 
found in section [VIII This main body of the text is supplemented with six appendices, where details concerning 
more computational or technical aspects have been gathered. 



II. GENERAL FRAMEWORK AND PRELIMINARY REMARKS 

The objects of our study are the scalar and vector form factors of the pion, defined through the following 
matrix elements 

(Tr t) (p 1 )7r a (p 2 )\m(uu + dd)(0)\n) = +Fjj° (a) 
(ir+(p + )Tr-(p_)\m(uu + dd)(0)\ri) = -F§(a), (II.l) 

with fh = (m u + rrid)/2, and 

i<7r + 7r|(B7 M u-d7 M d)(0)|n) = (p--p + )»F v (s), 

(II.2) 

respectively, in the presence of isospin-breaking corrections induced by the difference in the masses of the charged 
and neutral pions. In each case, s denotes the squared invariant mass of the dipion system, s = (p\ + P2) 2 or 
s = (p_ +p+) 2 , with Pi 2 — M 2 o, p\ — M%, and stands for the QCD vacuum state. The mass of the neutral 
pion is denoted by M w o , while M„ stands for the mass of the charged pion. We will define the isospin limit as 
the case when the neutral pion mass tends to the charged pion mass, M T o — > M n , while keeping the latter fixed. 
This explains the convention that we follow in this paper, namely that all quantities without superscript refer 
to the charged pion case (default case), and that we refer to quantities involving neutral pions by an explicit 
superscript. The minus sign in the definition of Fg(s) reflects a choice of phase for the charged-pion states. 
In addition, we choose the crossing phases to be —1 for charged pions and +1 for neutral pions. These choices 
are compatible with the Condon and Shortley phase convention in the isospin-symmetric situation. We further 
assume throughout that symmetry under charge conjugation holds. 

These form factors, while being perfectly well-defined observable quantities in QCD, are however not 
observables from a strictly experimental point of view: they can only be measured indirectly, and should thus 
at best be considered as pseudo-observables. For instance, in the Standard Model, the vector form factor Fy(s) 
appears in the physical process e + e~ — > h + ti^ through the exchange of a single neutral spin-one gauge boson, 
which in practice reduces to only photon exchange at low energies. As far as the scalar form-factors Fg (s) 
and Fg(s) are concerned, a similar statement can in principle also be made, but is of little use in practice, 
since the Standard Model contributions to the processes e + e _ — > ir + ir~ , 7r°7r°, arising from the exchange of 
a Higgs particle, are well below the level of sensitivity that one could expect for any experiment of this type 
in the foreseeable future. Despite these limitations on the experimental side, these form factors prove useful 
as a theoretical laboratory. They allow us to discuss and to illustrate several issues related to the structure 
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FIG. 1: Examples of radiative corrections usually included (left) and not included (right) in the data analyses. 

of isospin-breaking contributions within a rather simple context. The full complexity of experimentally more 
interesting situations, like the form factors or the decay amplitudes of light pseudoscalar mesons (eta or 
kaons) into three pions, can then be addressed on the basis of these considerations and the general framework 
developed here, see [22| and the forthcoming publication [3^ in the former or the latter case, respectively. 

In the present section, we aim at tackling two issues, namely discussing precisely the isospin contributions 
that we intend to deal with, and describing our general theoretical framework. Then we can focus on the specific 
pion form factors that we use as an illustration. 



A. Electromagnetic corrections 

At the fundamental level, isospin violations have two origins within the Standard Model: the electroweak 
interactions mediated by the gauge bosons, and the quark mass difference m u — m^, arising through the coupling 
of the light u and d quark flavours to the Higgs boson. Both effects contribute to the mass difference between 
charged and neutral pions. although the second one turns out to be marginal: the pion mass difference is mainly 



an electromagnetic effect [3J, [3 



Chiral perturbation theory lg, |36| including electromagnetism 37H42j| provides in principle a suitable 
framework to deal with these isospin-breaking contributions in the low-energy domain. It has bee n appli ed to the 
computation of several quantities at the one-loop level, including the 7T7T scattering amplitudes 



a appn 

mm 



43] and 



pion form factors [44( in the two-flavour case. Unfortunately, from a practical point of view, this is not a level 
of accuracy able to match the experimental one in several cases of interest (low-energy pion-pion scattering or 
K14 decay, for instance). One might of course contemplate the extension of this effective Lagrangian framework 
to next-to-next-to-leading order, but this is a more ambitious program, the interest of which might be limited 
eventually by the proliferation of low-energy constants. We will therefore not pursue this issue here. 



Instead, we will rather consider the point of view described in refs. 19. l2lj: we thus assume a situation, 
as is, actually, often the case in the analyses of experimental data (such as for E865 and NA48/2), where part of 
the radiative corrections due to real and virtual corrections have already been dealt with in some manner, while 
those that may remain are supposed to be negligible. In this kind of procedure, radiative corrections of the 
type shown on the left-hand side of fig. [1] for instance, together with emission of soft photons, are subtracted 
away, but other photonic effects, like the one shown on the right-hand side of the same figure, are not taken 
into account, but are considered to be negligible. Notice that the latter would be included in a genuine two- 
loop calculation in the framework of the QCD+QED effective theory. One might also think of considering the 
possibility of treating them in the dispersive framework that we are using here, for instance upon including also 
photons among the possible intermediate states in the unitarity conditions for the relevant partial waves. While 
this remains an interesting issue, it would however lead us beyond the purposes of the present work. Within the 
framework assumed here, this leaves the difference in the pion masses as the only remaining source of isospin 
breaking that we have to consider. In practice, it means that the charged and neutral pion masses will be kept 
at their experimental values, but the origin of the difference in their masses will not be addressed. In other 
words, we assume that general properties like analyticity, unitarity, and crossing, together with chiral counting, 
can be used to describe a world where the charged and neutral pion masses differ, even though the interaction 
at the origin of this difference is not explicitly accounted for. 
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B. Dispersive construction of the form factors 



The starting point of our study is provided by the dispersion relations satisfied by the pion form factors 
and scattering amplitudes. Here we will only be interested in the structure of these quantities in the low-energy 
region. In order to obtain dispersive representations that agree with their analytic structures up to two loops 
in the low-energy expansion, it is convenient to consider thrice subtracted dispersion relations (for a discussion 
of this issue in the isospin limit, see e.g. ref. 25| for the pion form factors, and ref. 23| for the 7T7t scattering 
amplitude). Let us start with the form factors, for which these dispersive representations read 



(0) 

mo) 
i 



l + l(r 2 r S "s + cfs> + uf( S ) 



F v (s) = l + ^(r 2 )y S + c v s 2 + U v ( S ). 



(II.3) 



In the last of these relations, the condition Fy(0) = 1, due to the conservation of the electromagnetic current 
and thus valid to all orders, has been used. Through crossing, Fg° (0) and (0) are equal to the corresponding 
sigma-term type form-factors, (ir°(p)\m(uu + dd)(0)\ir°(p)} and (tt^ (p)\fh(uu + dd)(0)\n :lz (p)), respectively, for 
which there also exist relations [45[ valid to all orders, that follow from the Feynman-Hellmann theorem 46 1. 



F§(0) 



m ■ 



dfh 



dMl 
dfh 



(II.4) 



Since the dominant contribution to the pion mass difference is purely of electromagnetic origin and independent 
of the quark masses 47[ , one has 



Fjm 

Ff(0) 



= 1 



(II.5) 



where the ellipsis denotes higher order terms. The unitarity parts are given in terms of dispersion integrals, 

- 3 i- dx lmFf(x)/Ff{0) 



U% (s) 
U v (s) 



TT 

t ,3 



x 3 x — s — iO 
dx lmFg(x)/Fg(0) 
x 3 x — s — iO 
dx lmFy(x) 
x 3 x — s — iO 



(11.6) 



In the low-energy region, the form factors are analytical functions in the complex s-plane, except for cut 
singularities on the positive real axis, starting at s = 4Af^ in the case of Fy(s), and at s = 4,M% in the cases of 
Fg(s) and Fg° (s). For s real and below these cuts, the form factors are real. In the chiral expansion, the form 
factors behave dominantly as 



ReF, 



00 - o(e% 



ImF, 



ir(ir°; 



(s)^O(E 4 ), 



ReF v (s) - O(E ), ImFy(s) ~ 0(E 2 ), 



(II.7) 



where E denotes either a pion momentum or a pion mass. Furthermore, intermediate states with more than 
two pions contribute only from the three-loop level onwards. Therefore, below the thresholds involving other 
states than the pions, and up to and including two loops in the two-flavour chiral expansion, only discontinuities 
arising from two-pion intermediate states need to be retained 25|, as illustrated in fig.[U In order to distinguish 
among the different tttt scattering channels, we use the following superscripts: 00 for 7r°7r° — > 7r°7r°, ++ for 



and tt + tt 



+ , H for 7T+7T 

7r 7T°. We have 



+0 for 7r + 7r — > 7T + 7r°, and x for the inelastic channels 7r°7r° 



7T 7T 



huFf '(-s) : = Re\lao( S )C(s)Ff*( S )e( S -4M 2 )-a( S )fS(s)Fr(s)9(s 



4Mi 



■0{E S 




FIG. 2: Diagrammatic representation of the unitarity relations for the form factors. In the case of the vector form factor, 
only the first diagram contributes. 



huF: (M = [-{<•<! a(s)ft(s)Fp(s)e(s - AM 2 ) - \ a (s)fZ(s)Ff*(s)6(s - AM 2 ) 



JmF^(s) = Re[a( s )fr(s)F™( s )9( s ~AM 2 )j+0(E 6 ), 
where we define the phase-space functions 

a (s) = Jl-^k, a(s) = 



0(E 8 ) 



1 



AM 2 



(II. 



(II.9) 



In these expressions, fo{s) and fi(s) denote the S and P partial waves, respectively, of the tttt scattering 
amplitudes A(s, t) in the corresponding channels. These partial waves have been defined as usual by projections 
of the corresponding amplitudes, 



1 f +1 
Ms) = — / dzA(s,t)P t 
S2ir J_ 1 



(11.10) 



with Pq(z) = 1 and Pi{z) = z the appropriate Legendre polynomials, and z = cos 9, where 9 denotes the 
scattering angle in the center-of-mass frame. The relation with the Mandelstam variables s, t, u (summing up 
to the squared masses of the incoming and outgoing particles) depends on the process under consideration. For 
processes involving four particles with the same mass, it is simply given by 

s - AM 2 



t 



(1-z). 



(11-11) 



This is, in particular, the case for A Q0 (s,t) (with M = M^o), as well as for A ++ (s,t) and for A + (s,t) (with 
now M — Mn). In the reactions involving both charged and neutral pions, it becomes 



~ (s - 2Ml a - 2Ml) + | V /( S -4M2 )( S -4M2) 



for A x (s, t), and 



2s 



(l-z), X(s) 



(M T + M,o) 2 ][ s -(M T -M w o) 2 



(11.12) 



(11.13) 



for A +0 (s,t). These expressions hold above the kinematical threshold, s > AM 2 for A x (s,t), and s > (M„ + 
M^o) 2 for A +0 (s,t), for instance. Let us point out that the normalization in pi. 101) differs from the usual 
definition of the 7T7T partial- wave amplitudes by a factor of 2: it would correspond to a decomposition of the 
scattering amplitudes given by 



A(s,t) = 16^y(2^+l)P,(cos0)/,(s), 



e>o 



(11.14) 



i.e. with the normalization factor 167r instead of the usual 32ir. This modification is motivated by the fact that, 
in the presence of isospin breaking, the two-pion states do no longer obey (generalized) Bose symmetry, except 
in the cases of two neutral pions, or of two identically charged pions. In these cases, the appropriate symmetry 
factor has been included instead in the expressions of the corresponding phase spaces, cf. eq. (jll.8l) . In the 
chiral expansion, the tttt partial waves behave as 



ReMs) ~ 0{E 2 ), lmh(s) ~ 0(E 4 ), £ = 0,1, 
Ref e (s) ~ 0(£ 4 ), lmf e (s) ~ 0(E 8 ), £ > 2. 



(11.15) 
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This is in perfect agreement with the chiral counting T[) of the form factors, together with the expressions 
pi.8[> of their low-energy discontinuities. 

C. Dispersive construction of tvtv scattering amplitudes 

It is actually possible to obtain dispersive representations for the two-loop 7T7t scattering amplitudes them- 



selves, following the same procedure as for the form factors. As in the isospin-symmetric case [23|, they follow 
from fixed-£ dispersion relations, combined with very general properties like relativistic invariance, unitarity, 
analyticity, crossing, and from the chiral counting properties that have just been recalled. Isospin breaking 
is not expected to modify the asymptotic high-energy behaviour of the amplitudes, so that two subtractions 
should be enough in order to obtain convergent dispersion relations. We start with three subtractions in order 
to construct low-energy expressions of the scattering amplitudes that are valid up to and including two loops 
in the chiral expansion. 

Whereas the various channels can all be described in terms of a single amplitude A(s\t,u) as long as 
isospin symmetry holds, several independent amplitudes, not related by crossing, are necessary in order to deal 
with all the different available channels once isospin is broken. Otherwise, the derivation proceeds as in the 



isospin-symmetric case 48(, up to the kinematical peculiarities due to the presence of particles with different 



masses. The relevant features can be inferred from the discussion of ref. 29(, devoted to the extension of the 
results of ref. Q to the scattering amplitudes of the mesons belonging to the octet of lightest pseudoscalar 
states, pions, kaons, and eta. We will therefore directly write down the resulting expressions, and then provide 
a few additional comments on their structure. 

i) Elastic scattering involving only neutral pions remains the simplest case, with a single, fully crossing 
invariant amplitude, which has the following two-loop structure (for convenience, we display, from now on, the 
dependence on the three Mandelstam variables s,t,u, although they are not independent) 

A 00 (s,t,u) = P 00 (s,t,u) + W°°(s) + W°°(t) + W™{u) + 0(E 8 ). (11.16) 

It involves a polynomial P 00 (s,t,u) of third order in s,t,u, symmetric under any permutation of its variables, 
together with a dispersive integral W Q0 (s). This function has a discontinuity on the positive real s-axis starting 
at s = 4:A'P , and specified by the I — partial- wave amplitude /q°(s), 

ImFy 00 (s) = 167rIm/ o 0O (s)6»(s-4Af2 o ). (11.17) 

Again, at the order under consideration, this discontinuity is provided by the unitarity condition, in terms of 
the £ = partial waves in the relevant channels, /o°(s) and /q (s), 

±-ImW™(s) = ^a (s) |/ 00 ( S )| 2 9(s - 4M^„) + a(s) \f£(s)\ 2 9(s - 4M^) + 0(E S ). (11.18) 

ii) The processes involving exactly two neutral pions, i.e. n^ir — > n^ir and tt + tt~ — > 7r°7r°, provide the 
next family of amplitudes that are related under crossing. They display the following structure at two loops in 
the chiral expansion: 

A x (s,t,u) = -P x {s,t,u)-W£(s)- [W+°(t) + 3(s-u)W+°(t)] - [W+\u) + 3(s - t)W+\u)} +0{E 8 ), 

(11.19) 

whereas A +0 (s, t, u) = — A x (t, s, u) through crossing. In the above expression, P x (s, t, u) represents a polynomial 
of third order in the Mandelstam variables, symmetric under exchange of t and u (Bose symmetry). The 
functions Wq\(s) and Wq(s) have discontinuities on the positive real s-axis, starting at s = (M^ + M n o) 2 and 
at s = 4M^ , respectively. These discontinuities are again given in terms of the appropriate lowest (5* and P) 
7T7T partial waves 

. 2 



ImW +0 (s) = 16tt 



Im/+ (s) + — * — '— Im/+ (s) 



](s-(M n + M n of) 



ImW 1 +0 (s) = 1671--^ Im/ 1 +0 (s)6»(s - (M n + M^f) 
X(s) 

Im^(s) = -167rIm/J(s)&(s-4M^o). (11.20) 
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Up to higher-order contributions, the unitarity condition allows to rewrite these expressions in terms of the same 
lowest partial waves. In the elastic channel, there is only one contribution, arising from the 7r + 7r° intermediate 
state, whereas the inelastic channel involves two contributions: 



^I^)| 2 + ^#IA + °( S )| 2 



1 ■lmW 1 +0 (s) = -i-- |/+°( S )| 2 6{s-(M n + M„ ) 2 ) + 0(E S ) 



l(s - {NU + M^f) + 0(E 8 ) 



16tt 1 w AV2( s ) 
JLhnWSW = -^ ( s )^( s )[/ ^ 

Hi) Finally, the subset of elastic scattering processes involving only charged pions remains to be considered. 
Their amplitudes being all related by crossing, it is enough to display explicitly one of them, for instance, 

A+-(s,t, u) = P + -{s,t, u) + [VK + -(s) + 3(* - u)W^-(s)} + [W+-{t) + 3(s - u)Wi~(t)] + W ++ (w) + 0{E S ), 

(11.22) 

where the third order polynomial P + ~ (s,t,u) is symmetric under exchange of s and t (Bose symmetry in the 
crossed u-channel). The three functions Wqi(s) and Wq + (s) have cut singularities along the real s axis, starting 
at s = 4M% or at s = 4M%. The corresponding discontinuities read 

lmW+-(s) = mnlmf+-(s)e(s-AM^ ) 
ImW?-(s) = l6Trlmf+-(s)9(s-4M%) 

lmW++(s) = 167rIm/+ + (s)(9(s-4Af 2 ). (11.23) 
The unitarity condition for the three tttt partial waves involved then leads to 

I i-ImW + -( S ) = a(s) \f+-(s)\ 2 e(s-4M*) + l* ( 8 ) \fS(s)\ 2 0{a - 4M W 2 ) + (D(E*) 
±.TmWr(s) = \fr{s)\ 2 e{s-AMl) + 0(E S ) 

-L W++( s ) = \a{s) |/ O + +( S )| 2 0( S -4M 2 ) + 0(E S ). (11.24) 

1D7T Z 

The reason why only the lowest S and P partial waves play a role in these expressions follows again from 
the chiral counting (|II.15[) for the partial waves. In the following, we will make use of the chiral expansion for 
the real parts of the I = 0, 1 partial waves, for values of s corresponding to the cut along the positive real axis, 

Refi(s) = <p t {a) + Ms) + 0(E 6 ), (11.25) 

with ip e (s) ~ 0{E 2 ) and ip e ( s ) ~ 0{E A ), so that 

\fe(s)\ 2 = [ReMs)} 2 + 0(E S ) = [^(s)] 2 + 2^(s)^(s) + 0(E 8 ), £ = 0,1. (11.26) 

Let us also emphasize that the functions W(s) only have a right-hand cut, that coincides with the right- 



hand cut of the corresponding S and P 7T7t partial- wave projections 24|. This structure is in agreement 
with the analyticity properties of the 7T7t scattering amplitudes A(s,t,u) required by unitarity and crossing. 
The decompositions (|II.16|) . (|II.19|) . and (|II.22|) satisfy these constraints, to the given order in the low-energy 
expansion. Of course, the partial-wave amplitudes have a more complicated analytical structure, coming from 
the projection in eq. 10|) . with also a left-hand cut, and even a circular cut in the case of the tt^-k — > tt^tt 
channel (for a description of the analytic structure of partial- wave amplitudes in a general context, see e.g. [3]). 
At this stage we should also stress that a full partial- wave decomposition (|II.14[) of the tttt amplitudes is actually 
not required. For our purposes, it is sufficient to know that the discontinuity of the latter in the complex s-plane 
can be written, in the low-energy region of interest here, as 

W(s,t) = 167r[Im/ (s) + 3^Im/i(s)] + ^> 2 (s,i), (11.27) 



with $£>2(s,f) ~ 0(E & ) as its dominant chiral behaviour. 
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over partial waves unitarity dispersion relation 



A at order E 2k 




/ at order E 2k 




Im / at order E 2h+2 




1 at order E 2k+2 
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4 — 











FIG. 3: Recursive construction of two-loop representations for the pion form factors and scattering amplitudes in the 
low-energy regime. One starts with polynomial expressions at order 0(E 2 ), k — 1, (resp. O(E ), k — 0) for the tttv 
amplitudes and scalar form factors (resp. for the vector form factor), and obtains the two- loop representations after two 
iterations. 

The very general features and the results that have just been presented are at the basis of the construction 
of two-loop representations of the pion form factors and scattering amplitudes in the low-energy regime. This 
construction is achieved trough a two-step recursive process of which we now give a short outline, summarised 
also in fig. [31 Chiral counting provides the initial information, namely that at lowest order the form factors 
reduce to real constants, to be identified with their values at s = 0, while the 7T7t scattering amplitudes consist of 
0(E 2 ) polynomials of at most first order in the Mandelstam variables. This initial input, together with unitarity, 
fixes the discontinuities of the form factors and of the amplitudes, through the expressions of the functions ipe(s) 
at next-to-leading order. The complete one-loop expressions are then recovered up to a subtraction polynomial 
of at most first order (second order) in s (in the Mandelstam variables) in the case of the form factors (of the 
scattering amplitudes). In turn, these one-loop expressions then provide the discontinuities at next-to- next- 
to-leading order, and thus the form factors (and amplitudes) themselves at order 0(E 6 ), up to a polynomial 
ambiguity of second order in s (third order in the Mandelstam variables). In the case of the 7T7t scattering 
amplitudes, crossing imposes further restrictions on the possible terms that may appear in these polynomials. 
Notice that the presence of these polynomials reflect the fact that the functions Wq°(s), Wq~i(s), etc. are 
only specified by their analytical properties, in particular as far as their discontinuities are concerned. This 
leaves room for polynomial ambiguities in the expressions of these functions, and the maximal degree of the 
polynomials is then limited by chiral power counting. 

This second iteration relies on the possibility to obtain analytical expressions for the 0(E A ) pieces ipt(s) 
of the real parts of the lowest partial waves from the one-loop 7T7t amplitudes, whose structures are no longer 
polynomial. This represents the technically most demanding step. The following sections are therefore devoted 
to the detailed implementation of this program. Since our main interest lies in discussing the effects of isospin 
breaking on the phases of the pion form factors, we will provide explicit two-loop expressions for the latter 
only. This means that we will stop in the middle of the second iteration of the recursive procedure. Completing 
this second iteration would provide full two-loop expressions for the form factors, and not only that of their 
imaginary part as needed for the phase shifts. This in turn would require one to obtain analytical expressions for 
the corresponding dispersion integrals, a daunting task as already explained in the introduction. We therefore 
defer this remaining step to future work. 

D. The subtraction polynomials: subthreshold parameters vs. scattering lengths 

The two-loop dispersive construction provides representations of the 7T7t scattering amplitudes that involve 
subtraction polynomials of at most third order in the Mandelstam variables. These polynomials depend on a 
certain number of parameters that are not fixed by the general properties, listed at the beginning of subsection 
III CI on which the representations for the irir scattering amplitudes rest. Furthermore, beyond general constraints 
coming, for instance, from the crossing property, there is, of course, a certain degree of arbitrariness in the form 
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of these polynomials, and thus on the physical meaning of the corresponding coefficients 



In the isospin-symmetric case treated in refs. 23l \24\, the form of the subtraction polynomial P(s\t,u 



was chosen so that some of its coefficients were identified with the subthreshold parameters of the amplitudes 
(the coefficients of its Taylor expansion around the center of the Mandelstam triangle). Among other reasons, 
this was motivated by the fact that the chiral expansions of these quantities show better convergence properties 
than, for instance, the scattering lengths. The latter can then be obtained from their expressions in terms of 
these subthreshold parameters (see, e.g., the corresponding two- loop expressions in Appendix B of ref. [24] and 
the discussion in ref. |17|). Along the same line of thought, in the presence of isospin breaking the most general 
subtraction polynomials of third order in the Mandelstam variables, and compatible with the symmetries of the 
amplitudes under crossing, can then be written as 



P°°(s,t,u) 



a 00 M\ 



F2 



3A 



(i) 



+ ^[(s- 2A4 2 o) 2 + (t - 2Ml a f + (u - 2Ml ?\ 

* TT 

o\(2) 

+ [ {s _ 2M ^f + ( t - 2Ml f + (u- 2M^o) 3 ] 

" TT 



+ ^-(s - 2Ml ){s - 2Ml) + %[{t- Ml Mt f + (u - Ml - M^f] 

TT TT 

A (3) 

+ %■[{*- 2Ml ){s - 2Mlf + (s 2Ml a f{s - 2M*)} 

" tt 

X (4) 

+ Jt[[t- Ml - Ml a f + (u- Ml Ml f] 



, v(2) 9i (2) 

+ ^ t( s " 2 ^") 2 + (* " 2M ") 2 ] + %(« - 2 ^) 2 
\(3) , x (4) 9X (4) 

+ ^6 t( S " 2M ') 3 + (* " 2M2 ) 3 ] + - 2 ^ 2 ) 3 - ( IL28 ) 

In each of these polynomials, the first line is of at most first order in s, t, and u, and corresponds to the tree-level 
amplitudes. The second line corresponds to the subtraction terms at one-loop level: to construct the scattering 



24 



amplitudes at one- loop precision, it is enough to consider only twice-subtracted dispersion relations [23 , 

On the other hand, it has become customary to rather let the scattering lengths play a prominent role, 
since they are usually considered to have a more direct physical interpretation than subthreshold parameters. 



The point we would like to stress here is that the framework developed in refs. 23, 24] , and that we extend in the 
present work to the isospin- violating situation, is rather flexible from this point of view, and can accommodate 
several choices of parameters, according to one's needs or purposes. It is simply a matter of appropriately 
choosing the forms of the lowest-order amplitudes and of the subtraction polynomials introduced at each of 
the two iterations. Different sets of parameters can be related, order by order in the chiral expansion, and 
the corresponding two-loop amplitudes differ only by higher-order terms. In the present article, we will use 
the representation in terms of subthreshold parameters, with the subtraction polynomials given in eq. f|II. 28[) 
above. In appendix [I] we provide the corresponding expressions in terms of the scattering lengths, and give an 
outline of how the present analysis can be reformulated in terms of the latter quantities. For more details, we 
refer the reader to our forthcoming work 
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Whatever choice is eventually considered, these polynomials altogether depend on fifteen independent 
subtraction constants. In the isospin limit, only six independent subtraction constants are required: isospin 
symmetry induces nine linear relations among these constants. These relations can be summarised by the 
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statements that, as M^o — > M n , one has 



where 



P 00 (s,t,u) 
P x (s,t,u) 



P(s\t,u) = A 



Ai_ 

F 4 



P(s\t, u) + P(t\s, u) + P(u\s, t) 
P(s\t,u) 

P(s\t,u) + P(t\s,u), 



2\2 



(s - 2M$ 



[{t-2Mlf + {u-2Mlf] 
[{t-2Mlf + (u-2Mlf] 



(11.29) 



(11.30) 



is the subtraction polynomial for the isospin-symmetric scattering amplitude A(s\t,u), cf. reference 24{. Some 
relations between these subtraction constants are given explicitly below [see, for instance, section UlI Al and the 
end of section IHI C| . 



III. FIRST ITERATION: ONE-LOOP EXPRESSIONS 



In this section, we discuss the pion form factors and the hit scattering amplitudes at leading order, and 
then proceed with the construction of the corresponding one-loop expressions along the lines described above. 



A. Leading-order form factors and mr amplitudes 

At lowest order in the chiral expansion, the form factors are constants, that may be identified with their 
values at s = 0, Ff(0), Ff(0), and F$(0) = 1. 

At the same order, the mr scattering amplitudes in the relevant channels read [cf. eq. 28[) ] 

A (M) = « " o M l - o M l° H 



Fl \ 3^3 71 J 3F2 
A + -(s,t) = !±-l s + t --MP > 



F% V 3 V 3F, 

A™(s,t) = (III.1) 

From these amplitudes, the partial wave projections are obtained as Kefi(s) — tpz(s) + 0(E i ), i = 0, 1, with 

2 , 2 , \ a x M 2 

-mi --Mir 





16ttF2 


Vo'( s ) = 


P+- ( 


32«F» { 




P+ - (s 


96ttF2 v 


= 


167T^2 ■ 



3 T 3 11 J 48ttF2 
-Ml 1 



3 V 247T.F2 
4M^) 



(III.2) 



The various parameters, like (3 X or aoo, that occur in these expressions are free, in the sense that they are not 
fixed by the general principles (analyticity, unitarity, crossing, and chiral symmetry). The occurrence, in these 
expressions, of Ml a rather than M% in the terms proportional to a x , a+_, or aoo is a pure matter of choice, 
and can be considered as part of the definition of these parameters. The presence, in the denominator, of F„, 
the pion decay constant in the isospin limit, is likewise a matter of convention. 
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As discussed in subsection III Dt in the isospin limit a x , cx f _, and aoo take a common value a. Similarly, 
the parameters f3 x and /3 + _ become equal to the same quantity j3 in this limit. Let us notice that there is 
no analogous quantity /3 o in the case of elastic 7r°7r° scattering, due to the Bose symmetry and the identity 



s + 1 + u = AMl . Both a and /3 were introduced in ref. [23| . They remain finite in the chiral limit, and describe 



the 1 = and 1 = 2 S-w&ve scattering lengths Oq in the isospin limit 23J, |24( at lowest order 



M 2 n M 2 

- (5a + 16/3), a 2 = _|- (« - 40). (III.3) 



l) 96ttF 2 y ^' ' u A8ttF 2 

On the other hand, the lowest-order S-w&ve scattering lengths corresponding to the amplitudes A x (s,t), 
A + ~(s,t), and A 00 (s,t) were computed in ref. [i^] and read 

2 / n 2\ / a a \ ~ M 2 

a 2 ) - (4/3 - a) 



3 v u u; v r ' 48ttF2 



,00 



As compared to ref. [40(, where (3 = 1 was taken at lowest order, we have kept the dependence with respect to (3 



in the isospin-violating correction terms. These expressions also account for the difference in the normalization 
of the partial- wave amplitudes as compared to [40 ] , see eq. (|II.14|) . Upon comparing these formulae with the 
expressions of the scattering lengths computed directly from the amplitudes displayed in (jIII.l|> . we obtain the 
following identifications: 

& = /3 + _ = 

M 2 -M 2 M 2 -M 2 
a x = a + 2/3 — 5 , o+_ = a + 4/3 — 5 , a o = a. (III.5) 

At higher orders, and in the absence of isospin symmetry, all these coefficients become independent, and these 
simple expressions receive additional contributions. The computation of the corresponding isospin-breaking 
corrections at next-to- leading order will be addressed below, see subsection IVI Al and appendix [El 



Finally, let us recall, from ref. 40], that the lowest-order irir amplitudes ([III. 1|) take the form 

A x (s,t) = -A(s\t,u) [s + t + u = 2M 2 +2M 2 ] 

A + -(s,t) = A(s\t,u) + A(t\s,u) [s + t + u = 4M 2 ] 

A ao (s,t) = A(s\t,u) + A(t\u,s) + A(u\s,t) [s + t + u = 4M 2 ], (III.6) 

with 

A(s\t,u) = S —^-, (III.7) 

and one has to be aware that the variable u that appears in A(s|t, u) takes a different meaning in each case, as 
indicated between brackets. Identifying these expressions with the ones in eq. then gives 

F 2 /. „ 2mB\ „ F 2 



at this order. Beyond leading order, the expressions (|III.8[) involve the low-energy constants 1% and £4 of ref. [lq . 
the appropriate formulae can be found in ref. [24|. At this point, it may be useful to make briefly contact with 
the discussion towards the end of subsection III Dt after eq. (III.28p . Indeed, one might actually consider three 
sets of independent quantities, (a, /3), (ap 1 , cLq), (£3, £4) as the unknowns of the problem. From a theoretical point 
of view, they are, to some extent, interchangeable. The last set naturally arises in the quark mass expansion that 
is implemented through the calculation of Feynman graphs generated by the effective chiral lagrangian. The 
two first sets are better suited for addressing phenomenological issues related to the analysis of experimental 
data. Here, we choose to organize the discussion in terms of the set (a,/3). The transcription in terms of the 
two S-wave scattering lengths can be found, as already mentioned, in appendix IFl and in a forthcoming article 
0. 
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B. Pion form factors at one loop 

We can now start the procedure described in figure [3] At this stage, the unitarity conditions take then 
the following form: 

lmFf(s) = \ a (s)^°(s)Ff(0)9(s - 4M* ) - a(s)<p x (s)F^(0)9(s - 4M*) + 0(E e ) 
ImF^(s) = a(s)<p+-(s)F£(0)9(s - - ~ a (s)<p* (s)Ff (0)9(s - 4M* ) + 0(E e ) 



ImF^(s) = ff(«K _ («)9(i-4M2) + 0(£ 4 ). 



(III.9) 



We have now to determine the full form factors, exploiting the fact that we know their analytic structure, i.e., a 
cut along the positive real axis, and the value of the discontinuity along this cut. We introduce the well-known 
functions Jq(s) and J(s) defined by the following dispersive integrals 



Ms) 

J(s) 



16tt 2 
s 



dx 



1 



4M 2 x x s 

dx 1 



iO 



a (x) 



167T 2 J4M 2 X X 



iO 



a{x). 



(111.10) 



These functions correspond to the standard one-loop integrals subtracted at ,s = 0, and through a change of 
variable the integrals can be brought into the more familiar form 



J(s) 



-1 



16tH 



dx In 



l-a:(l- a;) 



Ml 



(Hill) 



and a similar expression for Jo(s), with M n replaced by M^o. In the latter form, the integration is easy to 
perform for, say, s < 0. The expression of J(s) for the remaining values of s is obtained through analytic 
continuation, with the s + iO prescription on the cut, as made explicit in the representation (|IH.10[) . The result 
is well known and reads 



J(s) = 

Here, the function L(s) is defined by 



16tt 2 



[2 + <t(s) L(s) + ina(s)9(s - 4M%)] 



(111.12) 



L(s) 



< 




a(s)\ 


vl + 


a{s)) 










) + lJ 



[s > 4M 2 ] 
[s < 4M 2 ] 



1 - 



4Mi 



[s < or s > AMI] 



with er(s) = < 



4A /lM_i [0<s<4M 2 ] 



(111.13) 



where we have performed a similar analytical continuation of the phase-space function cr(s), whose cut extends 
over the interval < s < 4A1 2 .. The function ln(s) is defined as usual with its cut on the negative real axis. 
Analogous functions Lq(s) and cto(s) are defined upon replacing by M^o in the above expressions. Then 
one can easily find functions with the appropriate discontinuities pil.9j) to represent the form factors following 
the dispersive representation eq. (|II.3P : 



Uf(s) = P/(.s) + 16^i^ ( S )Jo( S )-167r^( S )-^^J» 



0(E 6 ) 



US(*) = Ps(s) + ^^-(S)J(S) - 167T 1 <f%(8) J (s) + 0(E 6 ) 

Uy(s) = P${s) + 16inpl-{s)J(s) + 0(E 6 ). 



(III. 14) 



Pg (s), Pg{s), and Py(s) represent calculable polynomials at most quadratic in s, which are determined by 
the property that the functions Ug v (s) and Ug° (s) have vanishing first and second derivatives at s = 0. These 
polynomials can be reabsorbed into the subtraction constants such as to build up the (one-loop) radii and 
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curvatures. At one loop, only one subtraction constant is required for each form factor. The corresponding 
expressions can therefore be rewritten as 



F§ (s) = Fg (0) 



F£(s) = F§(0) 



1 + a s s + 167r — - — J (s, 



16^(0) <p%(s)J(s) 



HmFS (0)-^(s)J ( S ) 



l + ags + 167npJ~ (s)J(s) 
F£(s) = 1 + + 16tt^-(s)J(s). 
At this order, the subtraction constants Og , a^, and a v are then related to the radii through 



„2\7r 



1 



Gal 



48tT 2 F% 

1 

96tt 2 .F 2 
1 



.F*°(0) 
Ff(0) 



Ml + M 2 M 2 
2/3, * - a, 



2f3 x 



Ml 
Ml + Ml 



Ml 



2 a °° 



4/3 + _ - 4c^_ 



Ml 



(r 2 ) v - „» v - ^^H- ■ 



while the curvatures are given by 



1 



2880tt 2 F 2 M 2 



5760tt 2 .F 2 M 2 



Ff(0) 
Ff(0) Ml 



M 2 Q 

28/3, - 2/3, — ^ + a 



28/3, - 2/3, 



M\ 

M 2 ' M 2 
A/ 2 



M 2 
00 IS 



M 2 



26/3+ _ + 4a + _ 



M 2 
Ml 



960tt 2 F 2 M 2 



(111.15) 



(111.16) 



(III. 17) 



C. One-loop representation of nn scattering amplitudes 



The form factors at two loops are obtained once we know the real parts of the one-loop S and P tttt 
partial wave projections, as shown in eq. (IIII.9I) . With this aim in mind, we now undertake the construction of 
the 7T7T scattering amplitudes to one-loop in the presence of isospin breaking. The starting point is provided by 
the lowest-order expressions (jIII.l[) of these amplitudes, supplemented with two amplitudes that are deduced 
from the former ones by crossing, and that are needed to express unitarity in the crossed channels, 



A +U (s,t) 



F 2 



t 



-M. 



J-2 



'-Ml, 



a x M 2 o 
3FI 



A ++ (s,t) = -% (s-^M; 



2c^-M 2 



Fl \ 3 V 3F 2 
together with the corresponding lowest-order partial wave projections, 



(111.18) 



16ttF 2 



A(f) 
2s 



-Mi 



-.Ml 



a x M 2 
48ttF 2 



Ae X(s) 



48ttF 2 2s 
16ttF 2 



-Mi 



a+-Ml 
24nFl 



(111.19) 



Applying the formulae given in sec. Ill C[ and recalling that at the one- loop order \fe(s)\ 2 = [^(s)] 2 + 0(E 6 ), one 
easily obtains the expressions for the unitarity parts of the various amplitudes, up to a polynomial ambiguity that 
can be reabsorbed into the corresponding subtraction polynomials P(s,t,u). For the amplitudes corresponding 
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to the elastic channels, with either only neutral ([11.181) or only charged pions (jll.241) . these expressions read 





= i[l6^°( S )] 2 J ( S ) + 


[16^§( S )] 2 J(s) 


w + -( s ) 


= [16vr^-( S )] 2 J( S ) + 1 


[167T^( S )] 2 J ( S ) 


W+-(s) 


= 36V S - 4M " )J(S) 






= \ [16^++( S )] 2 J». 





(111.20) 

For the amplitudes corresponding to the processes involving both charged and neutral pions (jll.20[) , one obtains 



Wg(s) 



(16^)V§(s)^ U (s)Jo( S ) - (167r)Vo + -W^(s)J(s) 



0i (M 2 -M 2 ) 2 



+ 



12F* 

PI 



9 

S Z S 



6(M 2 +Af 2 )] 
UM* - 14M 2 M 2 + 1LM* 



W+ U (s) 



3F 4 s 2 

/? 2 



,s--(Af 2 + Af 2 ) + 
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(M|-M 2 ) 2 



9Fi 



^+o(s) 



J+o( s ) 



36Fi 



2(M 2 + M 2 ) 



(M 2 - Ml 



2 \2 



J +0 (s). 



In these last expressions, we have introduced another one-loop integral [s + o = (M-k + Af^o) 2 ] 

'da; 1 X 1/2 (x) 



J+o(s) 



16tt< 



s+o 



x x — s — iO x 



together with the subtracted integral J+o(s) = J+o{s) — sJ^ o (0), i.e 



16tt' 



°°dx _ 

x 2 x — s — iO x 



1 



A 1 / 2 (a:) 



(111.21) 



(111.22) 



(111.23) 



The expression for J + o(s) can again be brought into the more familiar form of an integral over a Feynman 
parameter, 



JM = 



i 



16tt< 



dx In 



1 - 



x(l — x)s 



Af 2 - x(Af 2 - M 2 ) 



(111.24) 



through an appropriate change of variable. 

The functions Wq°(s), Wq~(s), etc. are defined by their discontinuities up to polynomial ambiguities. 
At one-loop order, these polynomials need only be of at most second order in the variables s, e, u. Taking into 
account the symmetry properties of the corresponding amplitudes, they may therefore be written as [see also 
eq. (|II.28|) and the discussion following it] 



P 00 (s,t,u) 
P x (s,t 7 u) 



+ 3A| [(s _ 2m1q)2 + {t _ 2Mlo)2 + {u _ 2Mlof] 



F 2 
A 



s - -M 2 - -M 2 
3 3 



(i) 



3F 2 

\(2) 



+ ^-(s- 2M 2 )( S - 2Af 2 ) + ^ [(* - Ml ~ Ml f + (« - A/ 2 - M 2 ) 21 



F 4 



F 2 



s+t- -Mi + 



xj 1 ! + A<. a i 



2A 



(2) 



[(* - 2Af 2 ) 2 + (t - 2Af 2 ) 2 ] + -^(« - 2Ml) 



(111.25) 
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The five subtraction constants Aqq , Ax , A+l, i = 1,2, are new free parameters. In the isospin limit, they are 
given by 

Ai + 2A 2 



AW 



A, 



A« 



A, 



A 



(i) 



(111.26) 



where, at this order, Ai,2 can be expressed in terms of the low-energy constants l\ and I2 of [16], cf. 24] and 
equation (|E.8j) below. 

The expressions for the irir scattering amplitudes that follow from the above results agree with those 
derived in refs. 



40. 
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431 ] from a Feynman graph calculation based on the low-energy effective lagrangian 



for QCD+QED, provided that one drops the contributions coming from virtual photons, while keeping the 
difference between charged and neutral pion masses [this procedure is discussed in greater detail in subsection 
IVI AI and in appendix IE] . 



D. Infrared behaviour of one- loop amplitudes and form factors 

In the isospin-symmetric situation, the hit scattering amplitude A(s\t,u), the vector form factor Fy(s) 
and the scalar form factor F£(s) are well behaved in the limit where the pion mass vanishes. For the pion 
form factors, this has been shown explicitly at the two loop level in jjjj]. This good behaviour in the chiral limit 
requires that some of the parameters that appear in these quantities develop themselves logarithmically singular 
terms in the limit whare the pion mass vanishes. This is necessary in order to compensate for the singularities 
coming from the unitarity part, for, as — > 0, 

1 / M 2 \ 1 
J ( s ) ~ 77TT ln — + FT- ( IIL2? ) 

'ifWO 16-7T 2 \ -S / 87T 2 

Taking into account the pion mass difference offers a wider range of possibilities from this point of view. Indeed, 
besides the path just described, reaching the isospin limit first, then letting the common pion mass vanish, two 
additional options might be considered, where one lets, say, the neutral pion mass tend to zero, keeping the 
charged pion mass fixed, or the other way around. In the first case, singular contributions in the unitarity part 
come from the function Jo(s), which behaves as in (|III.27I) . but with M n replaced by M^o. In the second case, 
eq. (|III.27[) applies directly. Notice that the function J + q(s) remains finite as either of these two limits is taken. 
It requires that both pion masses vanish for it to develop an infrared singular behaviour. 

Let us first consider the 7T7t scattering amplitudes obtained in the preceding sub-section. In order for the 
one-loop amplitudes A 00 , A x , and A + ~ to remain finite as M^a — > 0, with fixed, we must have 

a m M^ ->■ , a x M%„ ->■ a x M 2 , P x -> % , (111.28) 

and [for the sake of simplicity, we keep the notation for the pion decay constant, which has a regular 
behaviour in either limit] 

M 4 

a + _M%a * F2 a x (a x + 4l3 x ) In M* + finite 

M 2 - 

P+- ~ - 75 TB n Px{a x + ^ x )\nMl + finite 

A!, 1 ] -—j Pl In Mlo + finite , (111.29) 

whereas the remaining coefficients, Aqq'*, A+ 2 2, A^ 1 ^, and Ai 2 '* remain finite in this limit. Likewise, the coefficients 
a x and j3 x are free of infrared singularities. In order to avoid any possible confusion, we remind the reader that 
the coefficients aoo, a x , a + _ appear in the amplitudes multiplied by M^o as a pure matter of convention [see the 
remark after eq. (|III.lj) ]. Therefore, a x M^ and a + _M 2 need not vanish as M^o tends to zero with M v fixed. 
This feature is actually exhibited already by the lowest-order expressions given in eq. (|III.5j) . Furthermore, the 
quantities which appear on the right-hand sides of the above equations have to be understood as taking their 
lowest-order values. We have not distinguished them from their values at next-to-leading order, which appear 
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on the left-hand sides, in order not to overburden the notation. The appearance of infrared singular behaviours 
is a loop effect, and higher-order corrections will induce new singularities. At the next order, these may involve 
log-squared terms 25[, with an additional 1/(4ttF 7T ) 2 loop suppression factor. Concretely, from eq. 5 [) one 
obtains the lowest-order values 



a x = 2/3 , p x = (3. 



(111.30) 



Taking now the second limit, M n — > with M^o fixed, we find that the finiteness of the one-loop amplitudes 
A 00 , A x , and A + ~ makes the various coefficients [we denote their lowest-order, infrared finite, limiting values 
with a tilde on top of them, except for F^] behave as follows: 



A oo 

a x 

As 

A« 



48tt 2 F 2 
1 



48tt 2 
1 M* 



a x (a x +4# c )lnM* + finite 
/^lnM 2 + finite 

2a + _a x + /3 + _(S ; 



4a + _/3 K +/3 + _(a x + 4/y 



48tt 2 F 2 
1 M 2 
'96tt 2 F 2 

■-^P + _lnM' + finite 



hi Ml 
hi Ml 



finite 
finite 



5 M 2 



48?r 2 F 2 
1 Ml 



%-Pi_hiM< + finite 



^S + _/3 + _lnA/* + finitc 



A 



(i) 



n -8l hi Ml 
96tt 2 h+ - 



finitc 



A 



(2) 



1 



A^lnM 2 . + finite. 



48tt 2 

Now the lowest-order values inferred from eq. (|III.5|) read: 

a x = a -2(3, a + _ = a -4,0, fi x 



(111.31) 



(111.32) 



In appendix [E] we determine the various subtraction constants that appear in the expressions of the nir scatter- 
ing amplitudes obtained after the first iteration from the corresponding expressions obtained from a one-loop 
calculation. One may check that the formulae given in appendix [E] indeed exhibit the expected infrared be- 
haviour. 

Let us now turn towards the form factors, and consider the same two limits. As M w o vanishes while M r 
is kept fixed, the form factors remain free of infrared singularities provided 



1 



FJ (0) 



p x hiMl a + finite 



32tt 2 F 2 j?£(o) 
■ 7^2 ji %° (0) (S« ~ 2&) In Ml a + finite ; 



(111.33) 



while ag and Fg (0) remain finite [at this order]. Actually, in view of pil.30j> . this is also true for F£(0). In 
the case of the second limit, M„ — > and M^o fixed, infrared finite form factors require that 



1 



Fs (0) 



96tt 2 F 2 

1 F£(0) 
'16tt 2 F2 pf(0) 

1 Ml a 



(3 + _ hi Ml + finitc 

]} x In Ml + finitc 



48tt 2 F 2 



F£(0)(a x -2(3 x )lnMl + finite 
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a S ~ - 327r 2 F2 /3 + -lnM3 + finite 
1 M 2 ~ 

Fs(fl) ~ -^2-^ F 5(0)5 + -lnMj + finite. (111.34) 

Again, one may check that the explicit one-loop expressions given in appendix|E]reproduce the infrared behaviour 
obtained here from quite general arguments. 

To conclude this short discussion, we may observe that many of these infrared singularities disappear once 
the second mass also tends to zero, thus restoring the infrared features of the isospin-symmetric chiral limit. 
Upon studying the expressions (|III.16|) and (jIII.17j) in the light of the present discussion, a similar worsening 
of the infrared behaviour, as compared to the isospin-symmetric limit, can also be brought forward in the radii 
and curvatures of the scalar form factors. We thus can conclude that, to a certain extent, isospin symmetry 
tames the infrared behaviour of the soft pion clouds. 



IV. SECOND ITERATION: TWO-LOOP EXPRESSIONS 

So far, we have completed the first cycle of the procedure sketched in figure [3J The one- loop expressions 
of the form factors and scattering amplitudes obtained in the preceding section can now be used in order to 
construct the two-loop representations of the form factors. At next-to-leading order in the chiral expansion, the 
structure of the form factors is now as follows 



ReFS(s) = F§(0) [1 + TUs)} + 0(E 6 ) 



ReF£ (s) = Fg (0) 1 + T s (s) 



0(E 6 



(IV.l) 



and 



ReF v (s) = l + r£(») + 0(E A 



(IV.2) 



The one- loop corrections T s (s), Tg(s), and Iy(s) are easy to extract from the expressions of the form factors 
obtained in the preceding section, 



is(«) 



167r^° (s)ReJo(s) 



167r^+-(s)Re J(s) 



16^§( S ) 



Ff(0) 



Re J(s) 



167r-yg(s) F s f 0) Re J (s) 



T v {s) = a v s + 16irip+-(s)ReJ(s). 



(IV.3) 



As far as the discontinuities are concerned, they start at 0(E 4 ) for -Fg(s) and Fj$ (s), and at 0(E 2 ) 
for Fy(s). Using this power counting and eq. ([II. 251) . one obtains the discontinuities of the form factors at 
next-to-next-to-leading order 



ImFS (a) 



■ao(s)Ff(0){ l p° (s) [l + if (sj\ + iP° °( S )}9( S -4M 2 ) 
-<t( S )^(0)K( S )[1 + I s (s)} + r {s)}9{s-AMl) + 0(E 6 ) 



ImFZ(s) = a(s)F^0){^-( s )[l + T s (s)] + ^-(a)}0(a 
1 



4Mi 



' cr (S) 



)Ff(0) {^ ( S ) [l + lf(s)] + %(s)}6(s-4M* ) + 0(E e ) 



ImFy(s) = v(a){<pt-(8)[l + 11(a)] + tf-(8)}6(8-4MZ) + 0(E A ). 



(IV.4) 



A. Partial-wave projections from the one-loop amplitudes 



The computation of the one- loop corrections tpQ°(s), "0o( s )j V'o ( s )j i'l ( s ) t° ^he tttt S and P partial 
wave projections from the one-loop scattering amplitudes is less straightforward, and represents the next issue 
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to be addressed. In order to illustrate how this difficulty can be handled, let us start with the elastic scattering 
of neutral pions, i.e. the quantity ipQ°(s) describing the next-to-leading-order correction to the real part of the 
S'-wave projection for A 0Q (s, t) in the range s > 4M 2 . Trading the integration over the scattering angle for an 
integration over the variable t, with t_(s) = — (s — 4M^ ) < t < 0, we obtain (the functions J(t) and Jo(t) are 
real for t < 0) 

COO = ^^(Ss 2 - 16sM% + 20M*o) + ±- [I67r^°( s )] 2 Re J (s) + ±- [16tt^( s )] 2 Re J(s) 

1 1 (J2 A#4 pO II/" 

T " dtJ (t)+-—— T dt[16n^{t)] 2 J(t). (IV.5) 



167TS-4M2, F* J t _ (s) uw Stt s-4M* oJt _ {s) 

It turns out that the remaining integrals can be performed analytically. The relevant formulae can be found in 
appendix |A] The resulting expression can be written as 

16*<00 = 2^^—^^tf Q \ s )a (s) + 2t;g< \s)L (s) + 2^\s)^L v (s) 
+ 2^ ±) (s)a( S )ao( S )L( S ) + 2$ '°\s) (l ~ ^) L (s) 



+ ^°Hs)*o( S ) s M *° m2q Ll(s) + ^T\s)a (s) s _ A f M2o Z 2 v (s)j, (IV.6) 

with cr v (s) = a(s - 4M 2 + 4M 2 ) and L v (s) = L(s - 4M 2 + 4M 2 ). The various functions £$(s), £oo 0) ( s )> 
etc. that enter this expression of ipQ°(s) are polynomials in s and in the subthreshold parameters, which are 
given in appendix [Bj We have written the result in a way that allows for a straightforward connection with the 
similar expressions for the isospin-symmetric case, as displayed in ref. 2J]. Indeed in the limit M„o — ► M„ 
(and aoo, ot x , a + _ —> a, f3 x ,(3 + _ — > j3) one obtains the expected combination of I — and 1 — 2 contributions, 
weighted by the corresponding SU(2) Clebsch-Gordan coefficients, 



no, . . "oo, . _ o M ^ / ~ S 



F£ V s - 4M2 ^ 

" n—0 



k n (s), (IV.7) 



where, for the reader's convenience, we reproduce the expressions of the functions k n (s) of ref. |24| (the function 
&4(s) appears only in the P wave component ip^~(s), to be discussed below), 



1 / s - AM 2 1 

*»(«) = f 1 - — ) K.) • *3 (.) - 4- ; ■ 4= !»(.) , 
8vr V s J 16tt y/ s (s-4M%) 

For the remaining notation we refer the reader to [i^j (see in particular eqs. (3.36), (3.37) and the Appendix B 
therein). Let us, however, point out that the factor of 2 in (|IV.7|) takes care of the difference in normalization 
in the partial waves as compared to that reference, see eq. (|II.14[) and the comment preceding it. 

For the elastic scattering of charged pions, the computation of ipQ~(s) and of ipi (s), now in the range 
s > M%, proceeds along similar lines. The starting point is provided by the following formulae, 

i „ x (l) ,, x (2) 1 

= + - (--2A/J) 2 + A+ - 3 +f + - (s 2 -2sM„ + 4M*) + JL {16n^( S )f Re J (s) 



+i [i6^-( S )] 2 Rej» + ^j^J t _jm^m?m 



7T Jt-(s) 

dt Szz (* - 4M, 2 )(2 S + * - 4Mj) J(i), 



16tt S-4M2 7 t _ (>) 12Fi 
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and 



A?2 - A (2) 



1 



96^#" S(S " 4M ' } + 16. 36J* 
o 



' ^" ( S -4A/ 2 ) 2 ReJ(s) 



— n I dt\lQn^(t)} 2 ( 1 , 



2f 



+ 



1 



1 



32. s-4M2 Jt (s) 



{2[167r V J-(*)] a -[l^J + (t)] a } 



2f 



1(i,T .s- 41/_- ,/ t _( s) ^12^ 



rff • 



(£-4Af 2 )(2s + i-4M 2 ) 1 



s - 4Af2 
2t 



4A/ 2 



J(t) 

J(t), (IV.9) 



with now t_(s) = — (s — 4Af 2 ). Performing the remaining integrations with the help of the formulae given in 
appendix [S] leads then to 



M 4 



4A/ 2 



+<_fiw ( 1 - ^ ) l m + 2^( s ) ( x( S )ao( s )Lo( s ) 



+ 3£f;fj( s ) 



s 
M 2 



Vs(5-4M2) 



£ 2 ( S ) + 3d 3 -f2( S ) 



M 2 



^(s - 4M2) 



0" A (s) 
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16^r(s) = 2-f 



Fi V s- 4A/H 

(2;±), 



cr A (s) 



+ 2^( S )ll-^)l(,j 



s 
M 2 



v/ S ( S - 4M2) 
M 2 



i 2 ( S ) + 3^( S ) 



V S («-4M 2 ) 



1 M 2 
1 + — L(s) + fei 2 W 



(IV.ll) 



with er A (s) = <7o(s + 4Af 2 — 4M 2 ) and L A (s) = ^o(s + 4Af 2 — 4M 2 ). The various polynomials that enter the 
expression of ipQ~(s) and i/j^~(s) have been gathered in appendix [Bl Taking the isospin limit, as described above 
in the case of ipQ°(s), one recovers, for the S wave, the expected combination of / = and 1 = 2 contributions, 
and, for the P wave, the corresponding 1=1 component, 



j_ , , » , , , Mz I s 



E 

T n=Q 
4 



<T (-) - 2 ^ ^= ^ ^( S )M*), 

^ n— 



(IV.12) 



in full agreement with the results of 24 1 . 

Turning eventually towards the inelastic scattering ir + Tr~ — > ir°ir°, i.e. Vo( s )> the range of integration 
corresponding to — 1 < z = cosd < +1 is t_(s) < t < t + (s), with 



t ± (s) = -^{8- 2Af 2 - 2M 2 „) ± - yj( a - 4M 2 )( S - 4M 2 ). 



(IV.13) 



For s > 4Af 2 , one has t < and u < 0. In terms of an integration over t, one thus obtains 



\(!) 
Ax 

167rFi 



( S -2Af 2 )( S -2M 2 ) 



X (2) 

247r.Fi 



[ s 2_ s(M 2 +M 2 o)+4M 2 M 2 Q ] 
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FIG. 4: Examples of two-loop graphs of the factorizing (left) and of the non-factorizing (right) type involving pions with 
different masses. 

+ ^( S )^°( S ) ±- [2 + * (s)L (s)}+ ^"('KW " [2 + a(s)L(s)} 

ZTT TT 



1 
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8 ' f - V 


/{a - 4M2)( S - 4M2 ) ■ 


1 




24ttF4 


V /( S -4M2)( S -4M2 ) 
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/? 2 (M 2 -M 2 ) 4 


24ttF4 


v /( s -4M|)(s-4M2 ) 


1 





*+(«) 



dt 



4 I / - -A/r 



t+(s) dt 



Y (7i-6Af 2 -6M 2 ) 



3 



J+o(t) 
J+o(t) 



i+(a) 



f 2 



■>+(>(*) 



967rF - v /(s-4M2)(s-4Af2 o ) y*_ W 



/ dt{2s + t-2Ml -2M 2 ) 

(il^-M 2 ,) 2 



x i-2M 2 -2M 2 
With the help of the formulae displayed in appendix [5] one finds 

167^00 = 2 |f sjj^^ { &\s)a{a) + 2&Hs) 



J +0 (t). (IV. 14) 



1 



- 4Af 2 ) 

+ 2Ci 2;±) ( S ) (l - L(s) + 2^ 2 ^)( S )a( S )a ( S )i (.s) 

[£ 2 ( S )-£ 2 ( S )] I + 167rAi^(«) + 16ttA 2 ^(s 



A 1 /2 ( t_ (s ))£_ (s )_A 1 /2 ( i +(s))/:+(s) 



■3Ci 3) («) 



M 2 



^ S ( S - 4M 2 ) 



(IV.15) 



Here Ai^o(s) and A 2 V'o(' s ) represent two contributions that behave as 0(M 2 -M 2 ) and as O ((M 2 - M 2 ) 2 ), 
respectively, as one approaches the isospin limit. Their expressions, together with those of the remaining 
polynomials £i (s), etc. are given in appendix [Bl This expression of ipo(s) involves the function £±(s), defined 
as 



with 



xit) = 



£ ± (s) = \n[ X (t ± ( S ))}, 
y/(M v + M^o) 2 - t - y/(M„ - M^ f - t 

y/(M v + M w o)2 - t + y/(M n - M,o) 2 - t 



(IV.16) 



(IV. 17) 



Let us point out that the expression (|IV.15j) for ipo( s ) holds in the range s > 4M 2 , where the functions t ± (s) 
are real. An analytical continuation is necessary in order to describe V'o ( s ) m > sa yi the range 4M 2 < s < 4M 2 , 
as required, for instance, for ImFS (s), cf. eq. (IIV.4|) . For the applications that will be discussed in the following 
sections, we need not deal with this aspect, and the expression (|IV.15|) is sufficient. It is also useful to notice 
that in the isospin limit t+(s), A 1 / 2 (t + (s))£ + (s), and £|(s) all behave as O ((M 2 - M 2 ) 2 ). We keep however 
the contributions involving C + (s) as indicated in equation (IIV.15I) . so that each of the three pieces, when taken 
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separately, displays a regular behaviour as s approaches (from above) . Finally, in the limit where the value 



of the mass M_o tends to M,, one recovers the result of 



Mi 
pi 



3 

3 V 

s - mi ^ Q 



An) 
' ?2 



Jn) 
' SO 



fc n (s) 
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B. Two- loop representation of the form factors and scattering amplitudes 



Having obtained the partial wave projections in the S and P waves from the relevant one-loop tttt 
scattering amplitudes, we may now proceed towards obtaining the two-loop expressions of the form factors 
and scattering amplitudes. This requires one to evaluate the dispersive integrals in terms of which they are 
expressed. As mentioned previously, we will not be able to work out analytical expressions for all the integrals 
involved. Closed expressions will be obtained only for the contributions corresponding to so-called factorizing 
two- loop diagrams, see fig. |H They will involve the functions K n (s), defined as [24 1 



s r dx k n (s) 
K n {s) = - I — 



7T Jam" 2 *^ ' 



H) 



(IV.19) 



with the functions k n (s) given in eq. (|lV.8j) . and the understanding that Kq(s ) re mains denoted by J(s). 
Explicit expressions of the functions K n (s) in terms of J(s) can be found in ref. 24} . There will also appear 
functions K®(s), which arc defined in the same way in terms of functions k^(s), identical to the functions (IIV.8|i . 
but with the charged pion mass M n replaced by M n o. Similarly, we will keep the notation Jq(s) for Kq(s). 
Starting with the form factors, we obtain the two-loop representations 



F£ (s) = F§(0)ll + a%8 + b 



+ 8^(0)^°( S ) 



16^(0) <p%{ 3 ) 
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F i 



Ff(0)^T\s) - 2F£(0)^°)( S )1 [167T 2 J( S )J ( S ) - 2J{s) - 2J (s)} 



■A NF F£ (s) + 0(E 8 ), 
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and 



F£(s) = F^O^l + a^s + b^s 2 ) 



+ 167^(0)^" (a) 



1 + ClgS 



2 + , v 1 F S (°) x, s 
-ipt (s) 7-H- Vn(s) 



■8^(0)^(4 



tt Ff(0) 



J(s) 
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for the scalar form factors, whereas the vector form factor reads 



lGnip^ (s) 

Mi 



1 + a^s + ~<p$~(s) 



J(s) 



F 4 



7T 



+ A NP F£(s) + 0{E & ). 
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The contributions A NF Fg (s), A NF Fg(s), and A NF Fy(s) stemming from non-factorizing two-loop graphs, see 
fig. 21 are expressed as dispersive integrals, 
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A NF Fy (s) = 2 
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Their evaluation has to be performed numerically. Notice, however, that these representations are not neces- 
sarily best suited for this purpose, due to possible numerical stability problems. These functions are actually 
often expressed as two-dimensional integrals, which can be computed more efficiently. Examples of such repre- 
sentations can be found in the articles quoted under Q, and we refer the reader to them and to the papers 
quoted therein for further discussions on these aspects. 

We could proceed in a similar way in order to write down two-loop representations for the itit amplitudes, 
but this is not very useful for the applications of interest here. For the sake of illustration, let us consider the 
function W§ (s), involved in the amplitude for elastic scattering of two neutral pions, as an example. With the 
results already at our disposal, we obtain 

<°(*) = ~ [I6^g°( a )] 2 J (s) + [16^S( S )] 2 J(s) 
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The remaining W^-functions can be handled in a similar way, but we will not pursue the matter further here. 
As for the corresponding subtraction polynomials, they were already given in eq. (|II.28I) . 



V. ISOSPIN BREAKING IN PHASE SHIFTS 

In this section we turn to the issue of isospin breaking in the phases of the form factors, making use of 
the results obtained so far. We begin with a discussion of some general aspects of the phases of the form factors 
in the low-energy regime, and consider the lowest-order isospin-breaking corrections. The corrections at next 
order are then discussed in a framework where only contributions of first order in the difference of the pion 
masses are kept. 



A. General discussion and leading-order results 



The phases of the form factors are defined generically as 

F(s + iO) = e 2lS(s) F(s-iO), 



(V.l) 



where the phases will be denoted 5$ (s), Sq(s), and #f(s) for the form factors Fg (s), Fg(s), and Fy(s), 
respectively. For a discussion of the analyticity properties of the form factors, we refer the reader to ref. 21 1. 
Each of these phases 6t(s) has itself a low-energy expansion, Se(s) = 8ifl{s) + Se^s) + 0(E e ). Our aim is to 
address the issue of isospin-breaking corrections in the phases order by order in this expansion, i.e. our interest 
lies in the differences 



AS e (a) = 6 e {s)- &(«) = A 2 S e (s) + A 4 <5 £ (s) + 0{E 6 ) 



(V.2) 



between the phases Si(s) in the presence of isospin breaking and the phases 5i(s) in the isospin limit, with 

o 

ArA(s) = Si t n{s)— Si,n{s)- For the cases under consideration, we have 



Sf(s) = ±<ro(*)K°(*) + i,° °(s)]e(s-4M2 ) 
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■v(s) l * 



jg(0) 
Ff(0) 



(s - 4A/2) + 0(E 6 ) 



1 ,x*2(0) 



We thus deduce that 



i + rj(*) 

a(s) [rf-(s) + ^-{s))6{s-AMl) + 0(E e ) 



B{s-AMl ) + 0{E 6 ) 



Sl 2 (s) = <t(s)<pZ-(s)9(s-4A£) - ~a (s)<p x (s)8(s-4M* ) 
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FIG. 5: The lowest-order phase difference Aj 1 Sq (s) (solid line), defined in eq. ()V.9|I . as a function of the energy y/s in 
MeV, for a — 1.40 and /3 = 1.08. Also shown is the quantity A2S0 (s) (long-dashed line), defined in eq. (IV. 8j) . as well as 
the approximation to first order in A,r (short-dashed line), cf. eq. (IV. 151) . 



at leading order, while, at next-to-leading order, 



Ff(0) 



i + rs( s )-rs (s) 
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SoAs) = <j{s)^-{s)6{s-AMl 



~ 2 CT o( s )^o( s ) 
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(s-iM^o) 
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Here we have used the property that the quantities Fg° (Q)/Fg(0) — 1, rj(s), and Fg°(s) are all of order 0(E 2 ). 
At order 0(E n ), the isospin-symmetric phases Se, n (s) are obtained from the expressions (|V.4[) and (|V.5I) by 
setting MjO equal to M n , and by replacing the lowest order expressions of the 7T7t S and P waves, fi(s) and 

o o 

ipe(s), by their counterparts in the isospin limit, (Pt (s) and ipe respectively. In this limit, the differences 
Ff(0)/FJ(0) - 1 and rj(s) - Tf(s) vanish. 

Before proceeding with the actual calculation, let us make a few comments. First, we should point out an 
important aspect that emerges from the above expressions, and that has already been observed in ref. [2l|. At 
order 0(E 2 ), the phases of the form factors are entirely determined by the 7T7t scattering data. In the case of the 
vector form factor where, due to Bose symmetry, the 7r + 7r~ channel alone contributes, Watson's theorem is still 
operative: the phase of Fy(s) coincides with the phase of the P-wave projection of the corresponding scattering 
amplitude A + ~(s,t, u). For the scalar form factors, the situation is different, due to the mixing between the 
two channels that contribute to the unitarity sum in the S wave. Nevertheless, the phase has still a "universal" 
character, in the sense that its expression involves only the partial waves of the 7T7t scattering amplitudes, and 
makes no explicit reference to the form factors themselves. This property, however, rests entirely on the fact 
that Fg° {0)/Fg(0) = 1 at this order. Hence, this situation does no longer survive at the next order in the case 
of the scalar form factors. In addition to the universal parts, A^5q° (s) and A^Sq(s), provided by the O(E^) 
partial waves of the 7T7t scattering amplitudes, there now appear contributions AfSQ (s) and AfS^ (s) that 



20 



A 2 £o (mrad) 
10 r 




280 



300 



320 



340 



360 



380 



s ( MeV) 



FIG. 6: The lowest-order phase difference A\°Sq (s) (solid line), defined in eq. (|V.9|I . as a function of the energy yfs in 
MeV, for a = 1.40 and /3 = 1.08. Also shown is the quantity A2<5o (s) (long-dashed line), defined in eq. (IV. 8j) . as well as 
the approximation to first order in A,r (short-dashed line), cf. eq. (IV. 151) . 



depends explicitly on the form factors considered: 



A 4 <JJ (a) 



A u X7 ' 



X^Sq (s) + A* 8, 
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The universal parts A 4 <5q (s) and A 4 5J(s) correspond to the first lines of the expressions of 4 (s) and of 
<5q 4 (s) given in eq. (|V.5|i . respectively, from which the isospin-symmctric contributions are subtracted. The 
second lines in these same expressions correspond to the form-factor dependent contributions AfJJ (s) and 
Af 6q(s). Finally, we also note that for the scalar form factors some contributions in the expressions (|V.4[) 
and (|V.5|) start at s = M\, while others appear only for s > M%. This is of course the manifestation of the 
unitarity cusp in the phases themselves. 

At order 0(E 2 ), the isospin-symmetric phases 8e,2{ s ) are obtained from the expressions ()V.4[) . putting 
M n o equal to M„, and replacing the lowest order expressions of the 7T7t S and P waves y>i(s) by their counterparts 

o 

in the isospin limit, ft (s), which read 
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We then obtain 
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2(3 X -a x - Aa + _ A^ (f3 x -(3) + - /?) / _ 4 

96tt F 2 + 32nF2 ' S ' 



a-, - a Ml 
96tt F2 



a + _ - a Ml 



24tt F2 



cr( a )0( a - ml) 



A a aT(a) 
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0S+- - /?) (^-4M>( S ). 



If we make use of eq. (MI. 51) . these expressions become 



(V.8) 



A^sf(s) 



|^ h( S )0( S - 4Af2 ) - a (s)6(s - AMD] 
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32ttF2 
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and A% 6%(s) = 0. Here we have introduced the notation 

A, = Ml - Mlo 



(V.9) 



(V.10) 



Furthermore, Aj ^ (s) represents the leading-order isospin-breaking correction when the lowest-order relations 
flnLSj) for the subthreshold parameters aoo, ol x . . . have been used. However, the latter also receive higher-order 
corrections, that will be discussed below. This means that at next-to-leading order we also need to include a 
contribution 



Af°i/(s) = A 2 S e (s) 



aLO 
^2 



$i(s) 



(V.ll) 



which takes them into account. We show, on fig. [SJ a plot of A^Sq (s) as a function of energy [the numerical 
input values we use are given in table U and in equation (|VI.6[) ]. The curve exhibits the characteristic cusp-type 
behaviour 0,01 at the y/s — 2M^ threshold, where A2°<5J(s) takes its maximal value, close to 20 milliradians. 
For y/s greater than ~ 300 MeV, the value of A\°5q(s) stays practically constant, around 10 milliradians. 
The higher-order contributions mentioned in eq. (jV.llj) become sizeable only above the cusp. The curve for 
A\°8q°(s) follows a similar shape, see figure HI but its magnitude is reduced by roughly a factor of three as 
compared to A\°8q(s). In both cases, the difference A^ LO <5^(s) defined in equation (|V.11I) above is negligible 
in the cusp region, but becomes more and more important as the energy increases. 

It is also interesting to investigate the sensitivity of the isospin-breaking corrections with respect to the 
(unknown) parameters a and (3. This is done in figures |5] and [7] We notice that the dependence is strongest in 
the vicinity of the cusp, and is mainly driven by a. Indeed, at s = AMI, one nas 
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The value at the cusp of A^Sq" is directly proportional to a. In the case of A^Sq, the value is driven by 
the contribution proportional to /3, the contribution proportional to a being suppressed by a factor of three as 
compared to A\°8q . However, the relative variation in j3 only covers a restricted range, approximatively ±5%, 
so that the variations in a account for the largest part of the effect. As one leaves the cusp region towards 
larger values of s, the contribution proportional to a looses weight, and the variations are less ample. This 
behaviour is conveyed in a simple manner by the high-energy asymptotic expressions of A\°8q (s), and, even 
more strongly, of A\ 8^(s): 
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B. Isospin-breaking corrections at next-to- leading order 



The evaluation of Ai6^(s), the isospin-breaking effects in the phases at next-to-leading order, relies on 
the results obtained in the preceding sections. The corresponding numerical analysis will be the subject of 
section IVII below. Here, we wish to proceed for a while at the analytical level, but for simplicity, and since the 
next-to-leading order isospin-breaking corrections are expected to be small, we will restrict ourselves to the first 
order in A„. For this purpose, we expand the various quantities of interest with respect to A„, and neglect 
contributions beyond the linear terms. In the case of the phase-space factor for the neutral two-pion state, an 
expansion like: 



ct (s) = ct(s) 



f 



4M2 



A. + O(Al) 



(V.14) 



will not make sense when s remains close to AM 2 . This means that our expansion to first order in isospin 
breaking will only provide an adequate description in regions of phase space sufficiently away from the 7r°7r° 
and 7r + 7r~ thresholds. From an experimental point of view, this needs not constitute a serious drawback, since 
the vicinity of the two-pion thresholds is usually part of the regions of phase space where the acceptance is 
low, as can be seen, for instance, from 0, 0] in the case of the decay (see, however, the discussion on the 
— > 7r°7r°e ± i' e decay mode in [ijj]). From a practical point of view, we gain the advantage of having to 
deal with expressions which remain tractable. In the rest of this section, we will therefore remain within the 
framework set up by these two conditions - staying away from the two-pion thresholds, and considering only 
first-order isospin-violating effects. For illustration, let us consider the lowest-order corrections in the S'-wave 
phases, A£ <5J ( s ) and A£ <5J ( s ) that we have discussed in the preceding subsection. Applying the procedure 
that we have just described to the expressions (|V.9|) yields, for s > AM 2 , 
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On figures [5] and [6j respectively, these expressions for A£°(Jff(s) and A^ ^ («) are shown together with the 
exact expressions given in equations (|V.8|) and (|V.9|) . We see that the curve corresponding to the approximate 




FIG. 7: The lowest-order phase difference A 2 8q(s) for different values of a and f3. The solid curve corresponds to 
(a, 0) = (1.4, 1.08), while the cases (a, /3) = (1.0, 1.04) and (a, 0) = (1.8, 1.12) are represented by the long-dashed and 
short-dashed curves, respectively. 
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FIG. 8: The lowest-order phase difference A^ Sg (s) for different values of a and (3. The solid curve corresponds to 



(a,/3) = (1.4,1.08), while the cases (a,/3) 
short-dashed curves, respectively. 



(1.0, 1.04) and (a,/3) = (1.8, 1.12) are represented by the long-dashed and 



expression overshoots the exact formula by about 25% at y/s = 285 MeV, while the difference drops to a level 
around 5% for y/s above ~ 325 MeV. 

° 00 ° +- ° x 

The expressions of the next-to-leading partial waves V'o (s), ^o,i(s), and if>o(s) in the isospin limit have 
been given in eqs. (|IV.7|) . (IIV.12I) . and pV.18p . respectively. In order to evaluate A4<^(s), we will proceed as 
follows. First, we replace, in the polynomials Cqo\ s )i pi an d ^™'( s )j eacn occurrence of M% by M 2 . — A ff . 

Keeping only the terms of first order in A„ , we thus obtain 

= &\s) + ^§ S&H*) + O(Al) , (V.16) 

and so on. Next, we have to expand the functions that multiply the polynomials £,j^(s), Cl_.g, Cl_.p an d 

ix n \s), see eqs. (IIV.6I) . (|IV.10|) . (IIV.11I) . and (IIV.15I) . to first order in A ff . Finally, when subtracting from 

the functions £,oq\s), £+--Si C+l-p> Cx ^( s ) the corresponding combinations of polynomials £«*^(s) arising in 

the isospin limit, as given in (|IV.7j) . (|IV.12|) . and (|IV.18j) . isospin breaking only occurs through differences like 

a x — a, or /3 + _ — /?, for instance. Collecting these three contributions provides us with an expression for A4<5^(s) 

accurate at first order in A„. 

As an illustration, consider the case of ipQ°(s), cf. eq. (|IV.6|) . The first step is a straightforward algebraic 
— (n) 

exercise. It produces functions £ 00 (s), n = 0, 1, 2, 3, with 

CoO ( s ) = ^00 ( s ) + Coo ( s )> Coo ( s ) = Coo ( s ) + Coo ( s )> C ( s ) = CoO ( s ) + CoO ( V - 17 ) 

For the second step, we expand the various functions that appear in the expression for ipQ°(s) to first order in 
A v [the functions k n (s) are defined in equation (jIV.8|> ]: 

2^L ( S ) = 16^( S ) - 32vr^||fco( S ) + [4*o(«) + **(«)]} + O(A^) 

2 ^) Lv(s) = leTrfci^+STr^jfex^-fe^) - 7 i^ TJ [4ko(B) + hx(8)]\ + O(Al) 

2a(s)a Q (s)L {s) = 16^ 2 (s) + 8tt^ [k^s) - k 2 (s) ~ 4k {s)} + 0{A 2 n ) 

^^Ll(s) = IftrA^W - Wir^kw + ^i [3fc x ( S ) + 4fc 3 ( S )]] + 0(A») 
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3M>( S ) 



Ll(s) = 16nk 3 (s) + \^^-{- a [k l {s)-k2{s)\ 



Ml 



s-AM^ ~ vv " y " ^ Mr 

Collecting the two contributions allows us to rewrite eq. (jIV.61) as 

r4 3 
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[3fciOO + 4*a(*)U+ O(Al). (V.18) 
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k n (s) + 0(Al), 



where A£qo ( s ) ls the sum of <5£oo V s ) m eq. (|V.16[) and of the contribution generated by the expansions (|V.18[) 
to first order in A ff . The expressions of the functions AQq^(s) are displayed in appendix ID] Let us briefly 
comment on the appearance of contributions involving the factor M%/(s — AM%) in equation (|V.18|) . Upon 
closer inspection, one finds that the combinations [Ako(s) + fci(s)] /(s-4M|), [3fci(s) + Aks(s)] /(s-4M|), and 
— 4M 2 .), become actually proportional to a(s) as s approaches AM% from above. 
The extraction of the first-order isospin-breaking contributions from the remaining one-loop partial waves 
proceeds along similar lines, and we merely quote the resulting formulae: 
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More details, as well as expressions of the functions Ap_. s (s), A£|_.p(s), and A^x^(s) are given in appendix 
ID] Working at first order in A v has allowed us to cast the functions ipQ°(s), tpQi(s), and ipo( s ) m t° a form that 
makes the comparison with their expressions in the isospin limit straightforward. 

At next-to-leading order, the isospin-breaking contributions to the P-wave phase are thus simply propor- 

o _ 

tional to the difference i/j^~(s)— tp\ (s): 

4 r n , . a . . 1 

k n (s) + O(Al). (V.21) 
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In the case of the S'-wave phases, the corresponding corrections are naturally split into a universal contribution 
A4 5q(s) and a form-factor dependent piece Af #o(s), cf. eq. (IV.6I) . Keeping only the first-order isospin-breaking 
contributions, the two universal pieces read, again for s > AM%, 
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Concerning the form-factor dependent parts, the one-loop result (|IV.3I) gives, at the same level of accuracy, 
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VI. NUMERICAL EVALUATION 

This section is devoted to the numerical evaluation of the isospin breaking corrections A^Sg^s) keeping 
the full dependence on A ff . For this, we first need to know how the subtraction parameters that appear in the 
amplitudes and form factors after the first iteration are related to the corresponding ones in the isospin limit. 



A. Determination of the subtraction parameters 

From the dispersive representations of the form factors and scattering amplitudes, we have obtained the 
isospin-breaking corrections in the phases of the pion form factors beyond leading order. These expressions 
involve the normalizations -Fg(O) and F§ (0) and the two subtraction parameters a J and a" s in the one- loop 
expressions of the form factors, and only a subset of the 15 subtraction constants that appear in the irir 
amplitudes, namely aoo, ol x , a + _, /3 X , and on the one hand, Ago , A^ and A+_, i = 1,2, on the other hand. 
In the isospin limit, the latter are given, as indicated in equation (|III.26[) . in terms of the constants Aj discussed 
and evaluated in refs. 
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More accurate determinations have appeared since then [15j, ll7| , see below. 



(i) 

We thus merely need to evaluate the size of the isospin-breaking deviations like, say, Ai — A^. The subset 
aoo---/3+- is likewise related to the subthreshold parameters a and (3 in the isospin limit. At lowest order, these 
relations were given in eq. ()III.5|) . but in order to evaluate the phases at next-to- leading order, it is necessary 
to go beyond this approximation. Again, we only need to know the size of the deviations from the isospin-limit 
quantities a and j3. As discussed at the end of subsection IIII Al a and /3 represent the observables that we 
eventually would like to pin down from a phenomenological analysis of experimental data, so we have to trace 
down the dependence on these parameters beyond the lowest-order expressions. 

Let us now explain how we proceed with these tasks. For this purpose, we briefly come back to the 
discussion in subsection III Al The framework that we have presented there can be described by an "effective" 
lagrangian, whose form is similar to the chiral lagrangian used to treat electromagnetic corrections, but without 
including photons as dynamical degrees of freedom, as their effect is supposed to be treated by other means or 
otherwise to be negligible. The leading-order (strong) lagrangian £2 is then supplemented with a contribution 
of the form 



21]: 



£2 -> £2 + C{QUQU^) , Q = diag(2e/3, -e/3) , (VI.l) 

where C is a low-energy constant that breaks isospin symmetry among the pion masses. The last term, through 
its transformation properties under chiral symmetry, encodes the information about the electromagnetic origin 
of the pion mass difference. Although we could have absorbed it into the definition of C, we have left the electric 
charge e apparent, in order to make the comparison with the usual effective theory in presence of electromagnetic 
interactions more convenient. We call (|VI.1[) an "effective" lagrangian since there is no identifiable fundamental 
theory of which it would constitute the effective theory in the usual sense the quotation marks serving as 
a reminder of this limitation. Nevertheless, eq. (|VI. 1|) constitutes a suitable starting point for a low-energy 
expansion, with a well-defined and consistent power counting, which reproduces the features of the framework 
adopted here as far as isospin-violating corrections are concerned. Thus, the "effective" lagrangian at next- 
to- leading order is supplemented with the terms described in ref. [4o| . but with the corresponding low-energy 
constants denoted with a hat, to distinguish them from those obtained in the theory with virtual photons 
included. Indeed, the absence of virtual photons modifies the structure of the one-loop divergences, and the 
scale dependence of the renormalized low-energy constants hi- ((a) is given by 

e 2 ^kl{p) = -j^Zi, (VI.2) 
dfi lD7r z 

with (for the low-energy constants of interest in our case) 

1 A, . _ 1 A, 
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One has also a contribution quadratic in the difference M 2 — M 2 from the low-energy constant fc[ 4 (/i): 

We emphasize that these scale dependences are different from the ones of the equivalent low-energy counterterms 
discussed in ref. (4oT |. since we have considered a theory where no virtual photons arc included. Furthermore, 
they follow in a straightforward manner from the expressions given in eqs. (3.9)-(3.11) of [4pj j upon dropping 
the terms that do not contain the parameter Z = C/F 4 , with M% — M 2 = 2e 2 F 2 Z at lowest order, in the 
notation of that reference [in the case of CT14, only the terms in Z 2 must be retained, the constant k\^ being 
multiplied by e }. 

The relevant subtraction constants are then obtained upon matching the one-loop expressions obtained 
within the framework we have just described with the representations obtained in section lTTll The outcome of this 
exercise is displayed in appendix lEl Let us recall here that there exist explicit one-loop calculations of the various 



irir amplitudes 40j, |4l|, |43j and form factors [44| considered here, obtained within the full QCD+QED effective 



theory 37H42J. As mentioned in subsection III A| these calculations also include isospin- violating contributions 



arising from the exchanges of virtual low-energy photons, which are however not considered here. In order 
to make a comparison with these one-loop calculations, we must therefore remove the contributions of virtual 
photons from the expressions given in these references, and only keep the effects due to the difference of the pion 



masses 

,2 



2l| . From a practical point of view, this can be done as described above: contributions proportional to 



e , but without the appropriate Z factor, arise from the exchange of virtual photons and are discarded, while at 
the same time the low-energy constants are replaced by k\ (//). Furthermore, since we want to display the 

dependence on the two independent parameters a and /3, we have, in the computation of the loop contributions, 
explicitly kept the quantities F and rhB defining the leading-order amplitude (|III.7|) . for which we have then 
substituted the lowest-order expressions given in pil.8[) . This brings in another difference with the one- loop 
calculations available in the literature. 



B. Numerical input values 



We wish to investigate the size of the isospin-breaking corrections as functions of a and j3, for fixed values 



of Ai,2 and of the fc['s. As mentioned above, the former parameters have been evaluated before [15l Il7l 124 150| 
using sum rules and medium-energy tttt data. For the numerical evaluations below, we use the values from the 
"extended fit" of 



Ai = (-4.18 ±0.63) • 10 



-3 



A 2 = 



.96 ±0.12) • 10 



-3 



(VI.5) 



Let us notice that the sum rules that lead to this determinations of Ai and A2 also exhibit a mild dependence 
with respect to a and j3 24j,[50(. This dependence is, however, covered by the quoted uncertainties, and we will 
therefore not consider it further. 

As far as the constants k\ are concerned, we will assume that they take the same numerical values as 
the low-energy constants k\. Even though this identification constitutes an approximation whose precision is 
difficult to assess, it is not obvious to consider simple alternatives to this choice at the time being. Incidentally, 
this is also the option that was retained in ref. [19j ]. 

To obtain numerical estimates of the kjs, we proceed in several step. First, we make use of the relation 
between these two- flavour constants and their three- flavour counterparts K\ [37[, as worked out at one- loop 



level 



ref. 



52[. For the constants K\...K§ we then use the evaluation of ref. 53j, as given by the last 



line of Table 1 in this reference. These determinations rely on a set of sum-rules 54] involving QCD four-point 
functions, that are saturated by the lowest-mass resonances in the corresponding channels. This kind of minimal 



hadronic ansatz, which finds some justification in the limit of a large number of colors Nq [55l. |56|. is usually a 
good approximation j^]. We therefore endow the numbers of ref. 53] with a relative error of 33% (1/Nc for 
Nc = 3) accounting for neglected subleading effects in the 1/Nc expansion. We assign the same relative error 



to the constants K\ and K r xl that were estimated along similar lines in ref. 54[. Next, the relation between 



and also involves 



52j the SU (3) low-energy constants [57[ L\ and L r b . For the latter we have taken the 
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i 


1 


2 


3 


4 


5 


6 


7 


8 


14 


k\ (M p ) ■ 10 3 


8.4 ±2.8 


3.4 ± 1.2 


2.7 ±0.9 


1.4 ±0.5 


-0.8 ±6.3 


3.9 ±6.3 


3.7 ±6.3 


-1.3 ±2.5 


-0.4 ±6.3 



TABLE I: Values of the low-energy constants k\ used for the estimate of the subtraction constants. The values of the 
k\ correspond to the renormalized constants at the scale \i = M p = 770 MeV. The constants kz and kj do not depend 
on the renormalization scale, cf. eq. (|VI.3|l . 



0(p ) determination 10 • L\(M P ) — 1.46 ± 0.15 from ref. [58(. L\ is not so well determined, and can induce 
significant differences between the patterns of chiral symmetry breaking for two and three massless flavours 



5JJ-62j]. For our present purposes, we take 10 3 • L r 4 {M p ) = ± 0.5, a value which was advocated on the basis 
of the Zweig rule 57, 58 1, even though later fits and discussions favour larger values (sJSl. In any case, it 
turns out that h% plays only a minor role in the numerical evaluation of isospin breaking in the phases. Finally, 
the relation between the fc[s of interest to us and the -KJs also involve four constants, 8 9 15 for which there 
exist no reliable determinations. We have assigned an overall uncertainty of ±l/(167r 2 ) = ±6.3 • 10 -3 , based on 
naive dimensional analysis, to the values of 6 7 14 where these constants occur. The results of this analysis 
are displayed in table fl] Our values reproduce those given in ref. 52j , where however an overall uncertainty of 
±6.3 • 10 -3 was assigned uniformly to all the constants fc[(Af p ). 

For the sake of illustration, we let the parameters a and /3 vary within the intervals 1.0 < a < 1.8, 
1.04 < f3 < 1.12, suggested by the analysis of ref. 17J. These intervals cover a reasonable range of possibilities, 
but of course, if necessary, other values can be considered. 

To summarize, for the numerical analysis that follows we thus use as inputs the values of the constants 
Ai.2 in (|VI.5I) . the values of the constants fc[(yLt) as given in table Q] together with [fj^ 



Mir 



139.57 MeV, M_ 



134.98 MeV, 



92.2 MeV. 



(VI/ 



C. Size of isospin corrections to the phases 

Let us start with the phase of Fg(s). We will from now on restrict ourselves to the range of energies 
above the cusp. On figure [5] we show the isospin- violating correction ASq(s) = A 2 (5J(s) ± A4<5J(s), defined as 
explained after equation (|V.5|) . i.e., for s > 4M 2 , 



-a {s)ip x {s) 



F* s (0) 



-i - rs(s)-rs ( S ) 



(VI.7) 



The three terms of the decomposition in the second equality correspond, in succession, to the three lines of 
the first one, see the discussion after equation (|V.6|) . In the cusp region, A<5J(s) is rather well described by 
A 2 Sq(s) + A4 Sq (s), the contribution of Af Sq (s) is only marginal. At energies above ~ 300 MeV, Af Sq (s) starts 
to provide a sizeable negative contribution that more and more compensates for A^Sq (s), so that eventually 
A5q(s) ~ A 2 5J(s). This situation is reproduced if we take values of a and /3 different from the ones adopted 
for figure[51 but in the range considered here. We stress that it is A%Sq (s), and not A^ ^ (s), that provides a 
good description of the total effect in the region of higher energies. In this region A2<5J(s) and A2°<5J(s) differ 
by more than 2 milliradians, see figure [5] 

Turning next to the phase of Fg (s), the definition of the different relevant quantities reads 



^0 (s) 
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+ o*o( s ) C(*)-^o(«) C(«)-^oW 



Ff(0) 



= A 2 «JJ ( S ) + A 4 %^ (s) + Af «JJ («). 



(VI.8) 



We find a rather different situation than in the previous case. As can be seen on figure QUI the correction 
A4<5q°(s) is large, with A^Sq°(s) and A 4 5q°(s) both contributing in a substantial way. Here, A2<5J°(s), and 
even less so A^Sq" (s), do not provide a decent representation of the full isospin- violating contribution. Again, 
the situation shown on figureHUlfor specific values of a and j3 is actually generic for all values of these parameters 
in the ranges considered. 

Both A5%{s) and A<5J (s) receive contributions from iPq(s), which contains a piece A2"0q( s ) of order 
O(A^). We have checked that its numerical value is indeed tiny, in the range between —3 • 10~ 3 milliradian and 
—2 ■ 10~ 3 milliradian for all values of s and of the parameters a and j5 considered here. 

In the case of Fy(s), the form factor effects are absent from the isospin- violating contribution to the 
phase, which reads simply 



ASJ(s) 



A 2 5?(s) 



A 4 ^(s). 



a{s) (s)- 



fa- 



is) 



(VI.9) 



As can be inferred from figure the overall effect is small in this case, but the unitarity correction gives a 
substantial decrease of the lowest-order correction. 

In table UH we have also summarized our result in numerical form, for some values of the energy in the 
range between the 2M V threshold and the kaon mass. The presence of the form-factor dependent contributions 
in A4<5J(s) and in A^Sq (s) however precludes a direct application of our results to the experimentally more 
interesting case of the K^ 4 decay modes. This necessitates a dedicated study, which will be reported on elsewhere 
22| . Nevertheless, it is interesting to investigate the kind of conclusions that such a study might lead to from 
the analysis presented so far for the scalar and vector form factors of the pion. The quantity that comes closest 
to the observable of interest in the context of the decay mode — > 7r + 7r~e ± u e is the difference between the 
S and P phases, <5q (s) — Si(s), for which the total isospin-breaking correction reads 

A- t (s) = A6Z(s)-A5Z(s) 




FIG. 9: The total phase difference A5J(s) at next-to-leading order (solid line) for a — 1.40 and /3 — 1.08, as a function 
of energy (in MeV). Also shown are the lowest-order contribution A2#o ( s ) (lower dashed curve), and the combination 
A 2 5q (s) + A%S%(s) = AS^(s) - Af (s) (upper dashed curve). 
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FIG. 10: The total phase difference ASq (s) at next-to-leading order (solid line) for a — 1.40 and f3 = 1.08, as a function 
of energy (in MeV). Also shown are the lowest-order contribution A2<5q ( s ) (lower dashed curve), and the combination 
A 2 <5o°(s) + A%5f(s) = A5f(s) - A(8E°(s) (upper dashed curve). 



= A 2 ^( s ) + A^( s ) + Af<5 r.s)-A 2 ^( s )-A 4 <5 1 r ( s ). (VI.10) 

This correction in the phase difference is shown on figure 1121 together with the error band induced by the 
uncertainties on the various input parameters, for fixed values of a and /3. The main contributors to these error 
bars are the low-energy constants fe,, and in particular fci,2,34i which enter in the correction A4 Sq (s), through 
the isospin-breaking differences such as aoo — a, a x — a, and so on. Similar error bands have to be associated to 
the curves for A^S^s) and A^Sq (s) shown on figs. [9land[T0l respectively. Except in the vicinity of the cusp, 
the correction is thus relatively constant, around 11 milliradians for (a,/3) = (1.4,1.08), with an uncertainty 
that is somewhat less than ±1 milliradian. We show the correction, with the associated error band, for three 
sets of values for (a, (3). Despite the uncertainties, there remains a sensitivity with respect to these parameters. 

We have also compared the exact results obtained in this section with the approximation where only 
corrections of first order in A„ are retained, as discussed in subsection IV Bl For the range of parameters 
considered here, we find that using the approximate expressions for A^Sq (s), A^Sq (s), and A^J^s) does not 



(MeV) 


A 2 5S(s) 




A(5E(s) 


A8%(s) 


Aidfis) 


A%5f(s) 


A%5?(s) 


ASf(s) 


A 2 5?(s) 


A 4 S?(s) 


ASJ(s) 


286.07 


11.10 


1.42 


-0.29 


12.23 


3.53 


2.37 


0.38 


6.29 


0.007 


-0.001 


0.006 


295.97 


9.80 


1.24 


-0.38 


10.66 


2.70 


2.08 


0.61 


5.40 


0.027 


-0.005 


0.022 


304.89 


9.54 


1.19 


-0.45 


10.28 


2.41 


2.04 


0.78 


5.23 


0.050 


-0.009 


0.041 


313.47 


9.53 


1.17 


-0.51 


10.19 


2.27 


2.06 


0.91 


5.25 


0.077 


-0.015 


0.062 


322.02 


9.64 


1.15 


-0.57 


10.23 


2.21 


2.11 


1.04 


5.36 


0.107 


-0.022 


0.085 


330.78 


9.82 


1.13 


-0.62 


10.33 


2.20 


2.18 


1.16 


5.54 


0.141 


-0.030 


0.111 


340.17 


10.05 


1.11 


-0.68 


10.48 


2.23 


2.25 


1.28 


5.76 


0.180 


-0.041 


0.139 


350.92 


10.34 


1.06 


-0.74 


10.67 


2.29 


2.34 


1.41 


6.04 


0.228 


-0.055 


0.173 


364.52 


10.75 


1.00 


-0.82 


10.93 


2.41 


2.45 


1.57 


6.42 


0.295 


-0.076 


0.218 


389.71 


11.55 


0.80 


-0.95 


11.40 


2.70 


2.61 


1.85 


7.16 


0.430 


-0.126 


0.304 



TABLE II: The isospin-breaking corrections at order 0(E 4 ) to the phases (in milliradians) of the scalar and vector form 
factors for several values of the energy in the range between the 2M^ threshold and the kaon mass. The break-up into 
the various contributions discussed in the text is also given. 
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modify the values of A5q°(s), ASq(s), and A£J(s) by more than a few percents, as soon as the energy is 
more than ~ 20 MeV higher than the 2M W threshold, i.e.y/s > 300 MeV. For practical purposes, one possible 
option consists therefore in keeping the exact expressions for A 2 <5J°(s), A 2 Sq(s), and A 2 <5i(s), and in using the 
expressions truncated at first order in A,r for the next-to-leading contributions. 



A tot (mrad) 




FIG. 12: The isospin-breaking correction A^ ot as a function of energy, with the error band corresponding to the un- 
certainties on the input parameters, for (a, /3)=(1.8,1.04) (lower band, green), (1.4,1.08) (middle band, blue), (1.0,1.12) 
(upper band, brown). 
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VII. SUMMARY AND CONCLUSIONS 



In this article, we have addressed the issue of isospin breaking due to the difference between the charged 
and neutral pion masses in the pion form factors and ttt: scattering amplitudes in the low-energy domain. 
We have implemented a dispersive approach to obtain representations of the various tttt scattering amplitudes 
and pion form factors that are valid at next-to-next-to-leading order in the low-energy expansion. These 
representations rely on general properties such as relativistic invariance, analyticity, crossing, and unitarity, 
combined with the chiral counting for the form factors and partial waves, and provide an extension of the 
general frameworks developped previously in the isospin-symmetric case to the situation where isospin-breaking 
effects are taken into account. This construction needs as inputs the lowest-order expressions of the pion form 
factors and tttt S and P partial waves, and proceeds through a two-step iterative process, the partial wave 
projections obtained from the one-loop representation after the first step serving as inputs for the second step. 
At the two-loop level, we have obtained partially analytical expressions only, due to the difficulty of performing 
the dispersive integrals related to contributions of the non-factorizing type. We have nevertheless shown that 
in the limit where the pion mass difference vanishes, we reproduce known two-loop results for the scattering 
amplitudes and form factors in the isospin limit. On the other hand, we have obtained explicit analytical 
expressions for the phases of the two-loop form factors, on which we have focused. We have also provided 
somewhat more tractable expressions of the isospin-breaking corrections in the phases valid at first order in 
the difference M% — M% . This approximation provides a very good description of isospin-breaking effects at 
next-to-leading order in the phases, for all energies lying between ~ 300 MeV and the kaon mass. 

The dispersive representations involve a limited number of subtraction constants, which have to be fixed 
from experimental data or theoretical sources. We have related the subtraction constants involved in the phases 
of the two-loop form factors to their counterparts in the isospin limit, and we have provided a numerical 
evaluation of the isospin-breaking corrections in these subtraction constants. This has allowed us to perform a 
quantitative study of the size of isospin-violating effects in the phases of the form factors. 

We have displayed our results in terms of the subthreshold parameters a and f3 of the pion scattering 
amplitude in the isospin limit. These parameters represent the unknown quantities that have to be extracted 
from low-energy data. Equivalently, one may take the two S-w&ve scattering lengths Og and Og as unknowns, 
and we have provided the formulae necessary in order to operate the translation between the two formulations. 

As far as the phases of the form factors are concerned, the main difference with the isospin-symmetric 
situation is that Watson's theorem is no longer operative when the various tttt intermediate states that contribute 
to the unitarity sum become distinguishable, as a consequence of the explicit breaking of isospin symmetry 
through the pion mass-difference. The phase is still given by a universal contribution, expressed entirely in 
terms of data related to the tttt scattering amplitudes, but at next-to-leading order there appears a second 
contribution, that explicitly depends on the form factors under consideration. The numerical size of this form- 
factor dependent part is relatively small in the case of the tt + tt~ scalar form factor, but definitely more important 
in the 7r°7r° case. 

Wc have also investigated the sensitivity of the isospin-breaking correction to two subthreshold parameters 
a and /?, or equivalently to the two S'-wave scattering lengths Oq and Oq, which represent the quantities that 
have to be determined from data. Despite the uncertainties induced by the various other parameters, we find 
that the correction remains sensitive to the values of the subthreshold parameters. For the range of values that 
we have considered, and depending on the value of the energy, the effect in the total correction can represent 
up to 5 milliradians. The issue raised in the introduction about the possibility of a bias if isospin-breaking 
corrections are evaluated for fixed values of these scattering lengths remains therefore, in our opinion, open. 
Obviously, the situation that prevails after the present study devoted to the scalar form factors of the pions 
need not be representative of the one encountered in the case of the K14 form factors. To settle the issue, a 
dedicated study of this experimentally more interesting case is needed. This will be the subject of a forthcoming 
article [22 1. 
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Appendix A: Indefinite integrals of the scalar loop function 



In this appendix, we list the integrals involving the scalar loop function J(t) that are used for the 
computation of the one-loop partial wave projections in section IIVI For unequal masses mi =^ m 2 [when 
necessary, we assume that mi > 7712], the loop function J(t) is given by 



J(t) 



16tt 2 



In 



m 2 
mi 



+ ml 1 m 2 , A 1 / 2 ^ [t - A 1 / W _ (m 2 _ m 2 )5 



In 



77! 1 



2t 



In 



[t + AVa(*)]2 - (m\-m\f 



(A.l) 



where A (i) = t — 2t{m\+m2) + {m\— m|) ■ The above expression holds for t < (mi — m 2 ) . For mi = m 2 = M„, 
or mi = m 2 = M„o, this expression corresponds to the functions J(t) and Jo(t), respectively, defined in eqs. 
PILIOI) and Finally, the function J +0 (s) defined in eqs. (IIII.22|) and PIL24I) corresponds to the choice 

mi = Mrr and m 2 = M n o. 

The partial wave projections of the one-loop tttt amplitudes discussed require the computation of integrals 
of the type 



dtt n J(t) 



(A.2) 



where n can take negative or positive integer values. In the present context, the range —2 < n < 3 is particularly 
relevant. It proves convenient to define the variable x(tj) through 



t 



so that < % < 1 when — 00 < t < (mi 
replacing M T by mi, and M w o by m 2 , 



77! 7 



m 9 



77117712 X + 



1 



(A.3) 



7 2 ) 2 . The expression of x i n terms of t is given in eq. (|IV. 17|) . upon 



x(t) 



y/(mi + m 2 ) 2 - t - yj (mi - m 2 ) 2 - t 

y/ (777-1 + 777 2 ) 2 — t + \J (mi — m 2 ) 2 - t 



(A.4) 



For strictly positive values of the integer 77, the corresponding (indefinite) integrals take a relatively simple form, 
even when the masses mi and m 2 are not equal (irrelevant integration constants have been discarded), 



16tt 2 / dttJlt) = 



m 2 + m 2 m 2 
mf — mi mi 



(m 2 — m 2 ) In 



777 2 

mi 



i(t - mj - m 2 )A 1 / 2 (t) In X (t) + m 2 m 2 In 2 X (t) 



x t 



(A.5) 



16tt 2 / dtt 2 J(t) 



, 777 2 

In — 

i- 2 mi 



C 2 2m TO 2 

(7717 — m 2 ) In 

mi 



— (mi + m 2 + 10m 2 m 2 ) 



- [2t 2 - {m\ + m\)t - (m\ + m\ + lOm^)] A 1/2 (t) lnx(t) 



6 

+ m\m^{m\ + m 2 ) In 2 x(i) 



(A.6) 



16tt 2 / dtt 3 J(t) 



75 m 2 
m — 
7 2 mi 



(m 2 — m 2 ) In 



m 2 
mi 



12 
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— — (m\ + m 2 + STO^mf ) — — (mf + m 2 l )(m\ + to 2 + 28to 2 to,?,) 

- — [3t 3 - (m\ + ml)t 2 - (m\ + m\ + 8to 2 to 2 )£ 

-(m\ + m\)(m\ + m% + 28m?m|)] A 1/2 (i) lnx(t) 
+ m\m\[m\ + to 2 + 3to 2 to 2 ) In 2 x(t). 



(A.7) 



For n < 0, the corresponding expressions are more complicated, at least when the masses are different, but can 
be given a simpler form when expressed in terms of the function Hi t o(x) — — Li 2 (a;) — lnxln(l — x), which 
belongs to the family of harmonic polylogarithms 67]: 



I671- 2 / dtJ(t) = 



TO? + TO? TO2 

— 2 111 

m l — TO 2 TOl 

,2 2 



X t-X^t)ln X (t) + ^4^hi 2 X (t) 



+ (m( - m 2 ) 



Hi fi [ —X(t) ) - #1,0 I —X(t) ) + ln^i ln X (t) 



TO 2 



TO 2 



TOl 



TO-! 



TO 2 



(A.8) 



16tt' 
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AV2(<) 



2/ 2 



(m 2 — to?.) 2 



(m 2 — to 2 ) 3 

The range of integration, i_(s) < t < t + (s), with 
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(A.10) 



1 



t ± (s) = -- (s - 2ml - 2m 2 ) ± t ^ ( s - Am\)(s - 4to 2 ), 



(A.ll) 

will depend on the process under consideration. As the masses become equal, to 2 — } mi, one has [Ai 2 = ml — to 2 ] 



*_(«) = -(s-4to 2 )-2Ai 2 + 0(A 2 2 ) 



Then 



*+(*) = 

x-to = 

X+(s) = 

X(t_(s)) = 

X(t + (s)) = 



A 2 



0(Al 
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l-a(s) 
l + a(s) 
_ IA12 1 



(A.12) 



0(A 2 2 ) 



2 to 2 cr(s) 
s(s-4to 2 ) + Q(A 12 ) 
• A 2 2 + G(A 3 2 



Am' 



(A.13) 



where x±(s) = x(*±(s)). 



Appendix B: Polynomials of the next-to-leading-order 7T7t partial waves 

The expressions of the one- loop the partial- wave projections displayed in eqs. (|IV.6j) . (jiV.lOj) . (|IV.llj) . 
and (|IV. 15|) involve a certain number of polynomials whose expressions are given in this appendix. In the case 
of V>o°(s), these polynomials read 

= (5* 2 - 16M 2 oS + 20M* o ) 
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The polynomials involved in the expression for tp^ (s) read 
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while for ip^ (s) we obtain 
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-~/3 2 S 2 + ^ ( 5M 2 _ M 2 o ) _ ^2 ^g M 4 _ 6M 2 oM 2 _ 



1 



+ ^ x a x sAf 2 - p x a x Ml (2M 2 + M 2 ) 



1 



2304tt 2 M4 
1 

864tt 2 M4 
1 

864tt 2 M4 
1 

288tt 2 M4 
1 



4M 2 ) 2 



$l_M$ + ^ + _a + _MlMl - 2 a 2 _A^ + \pl_aMl 



288tt 2 M4 

Finally, for ipo( s ) the polynomials are given by 



Pi (-2M* + 2M 2 M 2 - 5 A^ ) + /3 x a x M% (2Af 2 - M 2 ) - M* 
[Pi {-AM* + 16M 2 M 2 - 25M* ) - 4/3 x a x M^ (2Af 2 - M 2 ) - a 2 x M* 



(B.3) 



A 



2A/4 

7T 

1 



(s - 2M 2 ) ( s _ 2M 2 ,) 



3MJ 



(s 2 - sAf 2 - sM 2 + 4A/ 2 M 2 ) 



864tt 2 A/4 

- s/3 x 



-/3 X (108/3+. 
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11/3 X ) + s/3 x A^ ( 36a + _ + — % + 27a 



91 



18/3+- (3M 2 + M 2 ) + — & (A/ 2 + M 2 ) 



9s/Va x M 2 



+ /3 X [p x (144M* + 144A^ - 112M 2 M 2 ) + 24£ + _M 2 (Af 2 + M 2 )] 

- 6/3 x A/ 2 (4a + _ + 3a 00 ) (Af 2 + Af 2 ) - 12/3 + _a x M 2 M 2 + 3a x M* (3a 00 + 4a + _ + 6a x )} 



1 



1 



2887r 2 Afi I 4 

,2 „ , r 4 



Pis 2 - P x s 



P x ( 3A/^ + -A^o 



+ a x A^ + /3 x a x Af 2 (2Af 2 - A/ 2 ) + /3 2 (7M* + lOM^b - 7M 2 M 2 „)} 



1 



128tt 2 M4 
1 

~ 1287r 2 Afi 
1 

8647r 2 Mi 



aooA/^o 



- |/9x (M 2 + M 2 ) + \a x Ml, 



+ ^Pl [I + -j^j (A/, 4 + M* - M 2 M 2 ) 
In addition, eq. (|IV.15[) involves two other contributions, 



P x s - 2 -P x {Ml + M 2 ) + ~a x M 2 



3/? 2 — (Af4 + A^„+Af 2 A/ 2 o) 



+ 2/3 x a x M 2 ( Ml - A/ 2 + ^ ) + a 2 M* ( 1 



s-- 3 m; 



Ml 



(B.4) 



16ttAi^(s) 



1 



967r 2 Fi 



/? 2 s + /3 x a x M 2 + /? 2 {Ml + Af 2 Q ) 



s-4M 2 



AV*(t_( a )K_( a ) 



and 
16ttA 2 ^( S ) 



1 



144vr 2 J F 1 i 



4M 2 ; 

, M 2 +A/ 2 0in A/, 



M 2 - M 2 M^o 



'^g + 1 W (s)) , +(s) 



(B.5) 



/3 2 S 2 - 5Pls{Ml+Ml ) 
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2/? 2 + 7M* - 6M 2 M 2 ) + 3/3 x a x M% s - 2/3 x a x M% (M% + M^ ) + 2a 2 x M* 
' ' ' [4$(5M* - 2M 2 M 2 + 5M*„) - 6/3 2 s(Af 2 + Af 2 ) 



+ ±p x a x Ml {Ml + Af 2 ) + 2a|M^ ] [J> (s) - («)] 
+ 6/3, [2/3 x (M 2 + Af 2 ) + a x Ml - [£ + (s) - £L(s)] 

+ 3/3 2 [H + ( S )-H_ (*)]}, 
where J-±{s) = J-(t ± (s)), and similar definitions for G±{s) and H±(s), with 

*±00 = 4(* - 2M| - 2M 2 ) ± i v /(.s-4M2)( s -4M2 ), 



(B.6) 



(B.7) 



and 

.F(f) = {Ml-Ml ) 

0(f) = (Ml + Ml )F{t) + (Mi-Ml )t\n 



In — — In yft) 



-(M 2 -M 2 ) 2 



(Ml + Ml, 2M n M„o X (t)) ^ - \ In 2 X (f) - ln X (f) + MtzMk m ^ZL 



4M 2 M 2 J r (i) - {Ml - M 2 )r In 



-(Ml-Mlo) 4 



(1 Ml + Ml M n \ 1 Ml - Ago M OT 
^2 M 2 - A/ 2 m M w o ; t J 2 M w c 



t 2 V(M2-M2 ) 



(B.8) 



Notice that individual terms in A2?/>q (s) may behave as C[A 2 xln(A Tr /Af 2 )] in the isospin limit, but cancellations 
occur between these terms, so that overall A 2 ^q{s) — 0(A 2 ). 

Appendix C: Two-loop form factors in the isospin limit 



In this appendix, we discuss the scalar and vector form factors at two loops in the isospin limit, where 
the complications due to the mass difference between neutral and charged pions are absent, and the dispersive 
integrals can all be expressed in terms of the known function J(s). The two form factors Fg° and Fg then 
become identical, since in both cases the two-pion states are projected on their / = 0, S-wave components, 
with identical Clebsch-Gordan coefficients. Furthermore, this exercise will provide a check of the calculation 
in the general case, which has to reduce to the expressions to be found below in the limit M^o — > M x . Let us 
recall that ref. 25| did not give analytical expressions for the form factors at two loops. Analytical two-loop 



expressions for the scalar and the vector form factors were given in [27] . The analytical expression for the vector 
form factor only had also been given earlier in ref. 26 1 . 

For the first iteration, the discontinuities of the form factors reduce to [in this appendix, we omit the 
superscript 7r most of the time] 



lmF s ( S ) = a( S )F s (0)cp ( S )8( S - 4Ml) + 0(E 6 ) 
lmF v (s) = a(s)ip 1 (s)6(s - 4Ml) + 0(E 4 ) , 



(C.I) 



with 



Ms) = <Po ( s ) " \ = ^o°0) - <Po(s) 



1 



167^ 



/3 ( s ~ - Ml \ + Z aM? 



tpi(s) = tp+ (s) 



1 

96ttF2 



P(s-AMl). 



(C.2) 
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The one-loop expressions of the form factors in the isospin limit read 

1 



F S (s) = F s (0) 



with 



F v 



U s (s) = 16tt<p (s)J(s) 



1 + ^)v» + 4 s 2 + U v (s), 
o 



Ml 



1 s 



16tt 2 F 2 36 Ml 



(8/3 - 5a) 



— (—)' 
360 \M*J 



(52/3 + 5a) 



U v (s) 



T^p{ImJ^ + ^ 2 ^ 



(C.3) 
(C.4) 

(C.5) 
(C.6) 



In order to implement the second iteration, it is necessary to include the next-to-leading contributions to 
the discontinuities of the form factors, 

ImF s ( S ) = a(s)F s (0){cp Q (s)[l + T s (s)] + M*)} 0(a - 4M 2 ) + 0(E 6 ) 
JmFv(s) = a(s){cp 1 (s)[l + T v (s)} + ^( S p( S -4M^) + 0(E 6 ). 



The relevant one-loop corrections Ts(s) and IV(s) to the real parts of the form factors are easy to obtain from 
their expressions given above, 



Ts(s) 



rv( s ) 



1 



/s 96^ 2 F 2 
1 



(8/3 - 5a) 



1 



<p (s)[2 + a(s)L(s)} 



(C.7) 



'v 



24tt 2 F 2 



+ -vji(«)[2 + <t(«)L(«)]. 



As for the one-loop contributions to the real parts of the S and P partial waves, they are conveniently 
expressed as 

Ms) = vrw - = \ - %{s) 

Ms) = M(s), 

where the expressions for ipQ°(s), tpQi(s), and tpQ (s) are given in eqs. (|IV.7I) . (IIV.12I) . and (|IV.18|) . respectively. 
They involve the polynomials ^\s) of reference 24 1. We reproduce them in a slightly different notation, in 
terms of the variable s instead of the relative momentum q = y/s/ 4M 2 — 1 in the center-of-mass frame, 



— i-j (105a 2 - 120a/3 + 392/3 2 ) + \ (HAi + 14A a ) 
4S2ir z 3 



1 864tt 2 



20 

(180a - 617/3)/3 - —(X 1 + A 2 ) } ^ 



311 „ /3 2 + - (UXi + 14A 2 ) 
1728tt 2 H 6 V ; 



^ a) (-) - 



1927T 2 

1 

11527T 2 



(a 2 + 4/3 2 ) 



5 -/3 2 ^ 



s 
7 



72tt 2 k M 2 576vr 2 



s 



(25a 2 - 80a/3 + 64^ 2 ) + 



1 



96tt- 



/3( 5 « " 8/3) 77? 



1 



M 2 32tt 2 



s 
Ml 



^ 3) ( S ) 



^ 0) («) 



5 (a 2 +4/3 2 ) + J-^/3 2 * 



288tt 2 
1 



864tt 2 



(93a 2 
1 



1728tt 2 



48tt 2 ' Ml 



8a/3 + 112/3 2 ) + -(Ai+4A 2 ) 
(207a + 256/3)/3 + I (X 1 + 7A 2 )| J- 
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+ {5H^ 2+ 5< A - +4A 4 G^ 2 

- sb<° 2 + ^ + 4 ' 32 » - sb«" + 2 « a?| + lab ^ 2 (s?) 

{i 0) (») = I ^ (15a 3 -460o,3 + 286^) 

+ {l72t ( "°- 12 ' 3) ' 3+ I< A '- A3 »}lJ| 

- -^« 5 °- 3 « + ™b« 5 « + «sj - n^ 2 (sj) 

1 ., (s-iMV ' 



2 



2 



72k 2 ' \ 4M2 

S! 4, (») = -jj^fa' + ^-Z^). 

Putting all elements together leads to 

1 3 
— -dmcFs(s) = Y,S n (s)k n (s)x6(s-4MZ) 

(C.8) 

4 

disc F v (s) = V n(s)k n (s) x 0(s - AMI) (C9) 



n=0 



with 



S n (s) = 16wip (s)6 nt0 
f 8tts 



<r2 >5 + 9iiif^- 5 °' 



+ 32ip (s) > <p (s)6 nfi 



+ 8[^o( S )] 2 ^ 2 + ||^ ) ( S ) 



V„(s) = l6Tripi(s)5 nfi 
f 8tts 
l~3~ 



+ 



< r2 >- + 24^2 13 



+ 32<^i (s) ^ <pi(s)5„ i0 



+ 8b 1 ( s )] 2 5„ ;2 + Ih^M. (c.io) 



Performing the dispersive integrals gives 

3 

U s (s) = Y,S n (s)K n (s) + P s (s) 



n=0 



(C.ll) 

4 

C/y( S ) = ^V n (s)tf„( S ) + Py(s). (C.12) 



n=0 



The two second-order polynomials Ps(s) and Py(s) are obtained as follows. First write 

2 3 r°° 



and, next, expand the polynomials 6>„(s) and V n (s), 



dx k n (x) 



x — s — iO ' 



Ml 
s 



VnW = vf + vw^ + vf) 



The polynomials Ps( s ) and -FV(s) are then given by 

n=0 
*■ n=0 



M 2 ) 



n=0 



Using the coefficients and k^ 2 ' 1 displayed in the following table: 
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(C.13) 



(C.14) 



n 





1 


2 


3 


4 




1 


1 


l 




1 


-4- 4 


1 


96 


16 


24 


192 


32 


576 T 


64 




1 


1 


7 




1 




19 


960 


192 


2880 


960 


128 


1920 ' 


3840 



we obtain 
Ps(s) 



P v (s) 



Mi 



16tt 2 F 2 
( Ml 



8/3 -5a ( s \ 52/3 + 5a / s \ 
36 \Ml) 360 \M%) 



\l6ir 2 F%J {\M%J 



ji^Sa 2 + + g!(a 2 + 0>) - ^(lUl + 14A 2 ) 



(135a 2 + 362/? 2 ) 



7r " - 2 „ fl2 , 87r 



M 



16tt 2 F 2 
M 2 



• — (a J - 2/? 2 ) + ^(89A! + 86A 2 ) 



9 \M^/ 60 \M%) 



(- 



Vi67r 2 F 2 y iv^ 2 y 

2 



-i-(45a 2 + 340a/? - 94/3 2 ) - ^(5a 2 + 50a/3 - 28/3 2 ) 
648 648 



Af 2 ) 



7 2 /„2\7T 



1 



6480 



(195a 2 + 1650a/? + 230/3 2 ) 



3240 



(10a 2 + 70a/? + 7/? 2 



16tt^ 



(Ai - A 2 ) 



The relations 



24] 



4M 2 \ - 1 

s _ 4Af 2 - - 3 - 

3L— r *K A (s)-3J(8)-±K 1 (8) 



1 



32tt 2 64 / M 2 



(C.15) 



(C.16) 
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allow us to eliminate K 2 {s) from Fs(s), and both K 2 (s) and Ks(s) from Fv(s), thus leading to the following 
expressions of the form factors 



F s ( S )/F s (0) = l + i(r 2 )5 S + c5 S 2 



+ 



+ 



Ml 



16tT 2 F 2 
2 



16tt 2 



\M 2 



P- -(8/3 -5a) )J(s) + 



1 s 
36 M 2 



{8/3 - 5a) 



360 ^^2; 



(52/3 + 5a) 



Vl67r 2 F2y I 9 



J(s) 



{w) ( 48 ^ 2 ( 11Al + 14A2 ) + 48^ 2 F 2 /3(r 2 )5 + i/3(551/3 - 15a)) 



+ 127T 2 ifl(s) 



S f ^ ) > - S f 1?3 ) - 3 «) + 4 (127/3 2 - 80a/3 + 10a 2 ) 



(^) ( 19m2 ( Al + A2 ) + 8 ^ F "(^ ~ 5a)(r 2 )5 



+ — (3684/3 2 - 1520a,3 + 25a 2 



+ 192tt 2 (11Ai + 14A 2 ) + -(976^ 2 - 480a^ + 285a 2 ) 



20 



27 \M 2 J 



27 \M 2 



27 v 



16?r2 r> , x 

+ —z-K 3 (8) 



(m 2 ) 



4 / A/f 2 

-27 (-f) (64/3 2 -80a/3 + 25a 2 ) 



/3 2 --(4/3 2 +a 2 ) 



{m 2 ) 



^(89Ar + 86A 2 ) + |^(8/3 - 5a)F 2 (r 2 )5 - ^(13322/3 2 + 135a 2 ) 



216 
16tt 2 



(2/3 2 - a 2 ) 



1 57T 2 

, (llAi + 14A 2 ) + ^(3224/3 2 - 2160a/3 - 45a 2 ) + ^(4/? 2 + a 2 ) 



-(8/3 -5a) 



■}■ 



(C.17) 



+ I le^j \ 9 



\mH {^ 2 ^ 2 Al ) - ^ + ^[167r 2 F 2 /3(r 2 )^ + 1] + 7 - 
(if) (t 48 ^( A2 - Al ) ~\] + \\^ 2 F 2 ^{r 2 y v + 1] + A/3(34/3 _ ll a) ) 
+ ^ (86/3 2 - 104a/3 + 9a 2 ) 



-4 



(if) 2 ^ 2 - 10 (if) m - a)+ 2(74/32 - 40a/3 + 5a2) - 128 (t" 

+ 87^4(5) (jj^j ^-^(^2) (44/3 2 - 10a / 3 + 5a 2 ) + ^(26/3 2 -40a/3 + 5a 2 ) 

o _ 



I P 2 



+ 



(m 2 ) 



1 ^f(\2 Ai) + |£/3F W V>£ + ^(H30/3 2 + 1650a/3 + 195a 2 ) 

-il^b~(71/3 2 + 35a/3 + 5a 2 
^-L(338/3 2 + 520a/3-45a 2 ) + ^ 



|^(52/3 2 -80a/3 + 10a 2 ) - ^/3 



(C.18) 
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In order to recover the expressions of [26| from these formulae, one simply needs to replace the various 
quantities by their expressions at leading or at next-to- leading order, as they can be found in refs. 16l. |24|. 

Ml 



r 



P 
Ai 
A 2 

„2\7T 
S 

2\7T 

IV 



1 



1 



48tt 2 
1 

48tt 2 
3 

1 

16tt 2 F 2 



32tt 2 F 2 
8^ 2 Fl [U 



(44 - 3£ 3 - 1) 



1) 



6 



12 
1) 



(C.19) 



At the end of this process, we then obtain perfect agreement with ref. [26 1. 



Appendix D: First-order isospin-breaking corrections to the one-loop partial waves 



The purpose of this appendix is to provide the explicit expressions of the functions that describe the 
isospin-breaking corrections to the one-loop partial waves, as given in eqs. (IV.19|) and (|V.20|) . 
In the case of tpoo( s ) the corrections that appear in eq. (|V.19|) read 

t| • A , ■ 2A,. ; ( " N ' 

1 Ml 



576?r 

j2 oQ„,fl i i i „ 2 



2 v 160^-+6^— 



2Aaf3 - 50A/3 2 - 203a ; 



72tt 2 s - AMI 



(AA(3 Z - 28a/3 + llcr) 



A^(s) 



1 



2304tt 2 

1 Ml 



12(32 w + 60/32 ^ + l2aP ip ~ 79a2 ~ 368/?2 



12tt 2 s - AMI 



1 



2304tt 2 
1 



' l2 P 2 T7l + 84 /3 2 T79 " 12a/3—— r - 53a 2 + 48a/3 - 64/3" 
Mi M 2 



(8/T - 12a/3 + 31a 2 ) + 



Ml 



(20/3 2 -4a^ + lla 2 ). 



(D.l) 



864tt 2 ^ ' ' 216tt 2 s - 4M 2 

In the case of tt + h~ scattering, we proceed as described in subsection IV Bl one first obtains the functions 
[X = S, P] 



>(i) / \ t(i;±) / \ . Ti 1 ^)/ \ t( 2 ) I \ ?( 2 ^) I \ , \ 

£+-;x( s ) = £+-;x( s ) + £+- ; x0) £+-;x( s ) = U-;x( s ) + f+-;xO) 
7(3) , s ? (3;±) , > ? (3;a) -(4) -(4;±) 7 (4;a). s 



(D.2) 



The expansion of the remaining functions gives 
a(s) 



2 

Mlo 



A w 



. L A («) - 16ttAi(8) - 8tt-^ [4fc (s) + fei(s) - fc 2 (s)] + G(A 2 ) 



V<s - 4M 2 ) 



Li(s) = 167rfc 3 ( S ) - Stt^ 



|fci(*)-|*a(«)+2*3(«) 



and 



M 2 



^ S ( S -4M 2 ) 



cr A (s) 



4M 2 



167r/c4(s) 
A 



-8tt 



M 2 



4/c 4 (s) - 



Ml 



AMI 



k 2 (s) 



(D.3) 



C(A 2 ). (D.4) 
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For the S and P partial- wave projections ipQ (s) and 1^1 this then leads to the expression (|V.20|) . with 
t(°) l o\ - 1 f_io«2 , 81 „2 _s_ _ 63 „ s „, R /1nfl 2_ 11(! 2 N 



A^ ;S (s) = 7^ f6/3 2 ^- 33/3 2 ^ + 15^-51a 2 + 16a/? + 56^ 



2304tt 2 V M£ ^ 2 ^ 2 



A ^ (S) = 2304^ ( v -6^ W + 39/?^-45a^-37a 2 + 64 a/? -48^ 

A £?- ; s00 = (26/3 2 - 2a[3 + 27a 2 ) 

A ^( S ) - 230^ ("1^ W + 8 ^ + ^- a2 - 4 ^-^ 

Agl s p(*) = T ^(26/3 2 - 10a/3 + lla 2 ) 

A ^ (s) = ^(5/^ 2 + 3a/3 + a 2 ). (D.5) 

Turning eventually to the inelastic 7r + 7r~ — > ir°ir° channel, one first rewrites the polynomials in eq. 
dTVUD as d n \s) = t n \s) + (AjFZ)5d n \s)+0(Al), with ^ 2) ( S ) = £? ;±) ( S ) + Next, one 



proceeds with the expansion of the remaining functions, 

a i/2 (m*)) r 1fl _ t 1C _ a. r t ,_ N , m 2 

^ - 4M 2 ) 



2 £_(«) = lforfciC*) - 16tt^ |*o(s) + ^jj^ [4*o(«) + *i(*)] \ + G(A 2 ) 



Af 2 A 
:£ 2 (s) = 16vrfc 3 (s) - 16tt 



AMI 



G(A 2 ). (D.6) 



Similar expressions with t_ (s) replaced by t + (s) are of order C(A 2 ), cf. equation (|A.13|) . and thus need not be 
retained in the present context. Finally, there are the two additional pieces to consider, 



-/3 2 -A? + a/3 + 2/3 2 
6 H Ml H H 



h(s) + O(Al), 



F* V s - AMI M l W 271 " 2 
A 2 ^ (s) = O(Al). (D.7) 

Putting the various parts together then leads to the expression given in (|V.20[) . with 

' l/^TTZ ~ 154 /3 2 T7? + 225a PlTF> + 352 ^ 2 - 270afi + 171a 2 



1728tt 2 \2 Mi ^ Ml r Ml 

1 Ml 



72tt 2 s - AMI 



(8^ 2 - 7 a/3 - a 2 ) 



A ^ 1)(S) = 23^ (-20/3 2 ^+3a/3^ + 80/? 2 + 36a/?+13a ; 

1 Ml 
+ 48tt 2 s - AMI 

A ^ 2)(S) = 23^ (-12/? 2 ^ + 57a/^ + 16/3 2 -104a/3 + 31a 2 

Is 1 A-f 2 

A ^ = 864^ W 2 Wl ~ ^ " ^ " 5 « 2 ) + 216^ T^lMl ~ <# ~ ^ 
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Let us close this appendix with a remark concerning the occurrence of contributions proportional to 
M 2 /{s — AM 2 ) in the expressions of the functions A£(")(s). When summed together into the functions ipQ°(s), 
ipo(s), and ipt~( s ) [they are absent in ipQ~(s)}, these singularities combine to give a regular behaviour as 
s — > 4:M%. Just like their lowest-order counterparts c/Jq°(s), </?§(s), ipQ~(s), and tp^~(s), the real parts of the 
partial- wave projections at next-to-leading order are regular at s = AM 2 , and the expansion in powers of A„. 
preserves this regularity, see also the remark following eq. (|V.19|) . 



Appendix E: Expressions of the subtraction constants 



The phases of the form factors discussed in section Ivl involve a certain number of subtraction constants, 
whose values are not fixed by the general properties underlying the dispersive relations that form the starting 
point of our construction. Two sets of parameters, aoo, ce x , and a + _ on the one hand, and /3 X and /3 + _ on the 
other hand, are directly related to the parameters a and /3 of the isospin-symmetric 7T7t amplitude, that are 
themselves related to the two scattering lengths in the S wave, cf. eq. 3[) . They represent the quantities 
to be extracted from experiment. What we need to know, however, is what becomes of the relations (|III.5|) at 
next-to-leading order. For the remaining set of parameters, the A's, the isospin-breaking corrections to their 
values in the isospin limit, given in eqs. (|III.26|) . also need to be worked out. In order to obtain this information, 
we have performed a one-loop calculation of the form factors and scattering amplitudes using the "effective" 
lagrangian approach described in subsection IVI Al The results of this calculation are shown in this appendix. 

First, at next-to-leading order, the expressions (|III.5j) become [the definition of the constants /C"° and 



K,2° in terms of low-energy constants introduced in 40j is given in eq. (|E.4j) below] 
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48tt 2 Fl 
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7/3(/3 + 5a) 



32^ 



(E.l) 



with = ln(M 2 / M 2 ) . Notice the occurrence of the term (/3 — a) -A?r i n the first expression. Since (3 — a ~ 
0(M 2 x InAf 2 ), this term reveals a logarithmic singularity (at most) in the chiral limit. Actually, it is finite as 
M^o — > 0. However, as M T — » 0, it develops an infrared singular behaviour, 



1 M 2 

^(7/3-4a)alnM 2 . 



32tt 2 Fl 



(E.2) 



We then obtain the following identification, at one-loop precision, with the various parameters involved in the 
polynomial part of these amplitudes: 
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From the one-loop expressions of the form factors, we obtain the following information on the subtraction 
constants in eq. (|III.15[) : at this level of accuracy, ag° and ay are unchanged as compared to the isospin limit, 
while 
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As far as comparison is possible, we find agreement with the existing results in the literature quoted at the 
beg inning of this appendix, except in two instances. The expressions for charged pion scattering given in ref. 
43j only included corrections of first order in isospin breaking, with which we agree. The formulae we give here 



51 



are not limited to this approximation. Furthermore, we found a slight disagreement with the result of 4J| for 
Fy(s): the radius (r) v exhibits an infrared divergence proportional to lnAf 2 as M^o goes to zero, whereas 
we find that (r) v remains finite in this limit, but diverges as ~ lnAf 2 if we send the charged pion mass to 
zero, keeping M n o fixed. This is also what follows from our analysis in subsection IIIIDI [69]. In this context, 
it is important to stress that in the expressions given above, the scale-independent low-energy constants ti are 
defined as 



i.e. the normalization of the logarithm is provided by the charged pion mass. We have also checked that the 
results given in this appendix display infrared behaviours in agreement with the ones obtained in section fill Dl 
provided one takes 

Ai = i(V^V, A 2 = -^(Vf)/? 2 - (E.8) 
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Notice that these expressions differ from the ones given at the end of appendix [Cj see eq. (|C.19|) , by the factor 
ft 2 . Both are compatible at one-loop order, where one would take /3 — 1 in the above formula. 



Appendix F: Two-loop phases in terms of scattering lengths 



It is perfectly possible, within the framework adopted in this article, to write down expressions that involve 
the scattering lengths instead of the subthreshold parameters. This is achieved by choosing a parameterization 
of the lowest-order amplitudes in terms of the scattering lengths, i.e. the value of the amplitudes at their 
respective thresholds, rather than in terms of their values at the center of the Mandelstam triangle, as done in 
the rest of the present article. The expressions 1 [) and (|III.18I) are thus replaced by 
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At this order, the relation between the two sets of parameters is simple, 
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The quantities a — a x ,a + ^,aoo etc., are scattering lengths to the extent that the tree-level amplitudes (jF.lj) 
satisfy 



Re^4(s, t, u) | h = Wira. 



(F.3) 
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The parameters a will keep their meaning up to next-to-next-to-leading order if the above relation still holds 
for the two-loop amplitudes. This can be achieved upon adjusting the subtraction polynomials accordingly. In 
practice, this is done through the following choice: 
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where we have subtracted the values of the one-loop integrals at the appropriate kinematical points to ensure 
eq. (jEU) 
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These expressions should involve the same number (fifteen) of independent subtraction constants (among them 
now the scattering lengths) as the ones given in eqs. (|II.28[) . This means that there exist six relations between the 
twenty-one parameters occurring in the above polynomials, which stem from crossing symmetry [by construction, 
the unitarity parts satisfy separately the crossing relations], P x (t,s,u) + P +Q (s,t 7 u) — and P + ~(u,t,s) — 
P ++ (s,t,u) = 0. This yields 
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and 
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For the time being, it is convenient not to make use of the two last relations, and to treat all the ,9- wave 
scattering lengths as independent. Notice also that in the chiral counting the scattering lengths are of order 
0(E 2 ), whereas b x , b + _ and b ++ are of order O(E ). One may now repeat the computation of the relevant 
partial-wave projections starting with the expressions of the lowest-order amplitudes in terms of the scattering 
lengths a,i and effective range parameters 6, (i = 00, ±0, x, + — , ++), considered as independent quantities, 
following the procedure outline in sec. [TT| and fig. [3] The results can still be brought into the representations 
(|IV.6j) . (|IV.10j) . (|IV.11|) . or (|IV.15j) . but the expressions of the polynomials involved are different from the ones 
given in appendix [B] For the scattering of neutral pions, the polynomials for ipQ°(s) now read 
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For the scattering of charged pions, the polynomials involved in the expression for ipQ (s) (S'-wave) read 
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while for the P-wave contribution, tpf (s), we obtain 
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Finally, in the case of a scattering involving two neutral pions and two charged pions, we obtain the 
expression for the polynomials describing tpQ 
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In addition, eq. (|IV. 15|) involves two other contributions, one of order 0{A 7T ) which reads 
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We do not give the explicit expression of A2'0o, since it represents a tiny contribution of order 0(A 2 ) which 
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can be neglected for practical purposes, as indicated in section fVTl 

At this stage, one can follow the discussion of section [V] and determine the isospin-breaking differ- 
ences ASq, ASq° , A$i in terms of the different scattering lengths and effective range parameters. Even 
though one might hope to determine all these parameters from high-precision data on the different chan- 
nels involved, it seems more realistic to express them in terms of the subthreshold parameters oti, ft, \^ with 
i = 00, ±0,2;,+-,++ 
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These expressions can be exploited, by relying on appendix [E] and expressing the subthreshold parameters 
ctfj, pi, \\ with i = 00, ±0, x, H — , ++ in terms of the isospin-limit parameters a, /3, A^™^ . The latter could be 
taken as the fundamental parameters of the analysis, but they can also be traded for the two 7T7T scattering 
lengths eig and a 2 , (up to higher-order corrections that can be estimated using Chiral Perturbation Theory). 
This series of matching will be indeed the point of view adopted for the analysis of decays, allowing us to 
reexpress the isospin-breaking correction to be applied to the phase-shift difference in terms of the two scattering 
lengths Oq and a 2 , [22|. 
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